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Zusammenfassung —EKs wurde die BlasengrOsse in Blasensiiulen an einigen Modellsystemen 
untersucht. Dazu wurden die Blasenschwiirme in verschiedenen Hohen oberhalb der Gasverteiler- 
platten fotogratiert. Die Hiufickeitskurven fiir die ideellen Blasendurchmesser stellen in erster 
Niherung svmmetrische, eingipflige Verteilungskurven dar, die sich durch einen mittleren 
Durchmesser d, und die Standardabweiehung o charakterisieren lassen. Die Filterporengrésse 
der Verteilerplatten beeinflusst die Grosse der entstehenden Blasen nur sehr wenig. Mit wachsen- 
dem Durchsatz steigt an einem festen Ort in dev Saule der mittlere Blasendurchmesser q., bis 
zu einem Grenzwert an. Mit wachsender BeobachtungshOhe in der Siule nimmt unter sonst 
gleichen Bedingungen der mittlere Blasendurchmesser d,, zu. Von verhiltnismiissig geringem 
Einfluss auf die Blasenyrésse ist der Siulendurchmesser. Das Zusammenlaufen der Blasen beim 
Aufsteiven in der Saule lisst sich durch Elektrolytzusiitze oder besser durch geringe Zusiitze 


obertlichenaktiver Substanzen weitgvehend vermindern 


Abstract Bubble sizes were investigated in bubble columns with the aid of some model systems. 
Bubble swarms were photographed at various heights above the gas distributor plates. The 
frequeney curves for the ideal bubble diameter show to a first approximation a symmetrical, 
single peaked distribution characterized by a mean diameter, = and a standard deviation ce. 
The pore size of the distributor plates has only a very slight influence on the size of the bubbles 
produced, As the yvas-flow increases the mean bubble diameter, Goan increases up to a limiting 
value at a fixed point in the column. As the observation level is shifted upwards the mean bubble 
diameter, d,,. increases, other factors remaming constant. The diameter of the column has a 


relatively slight influence on bubble size. The coalescence of bubbles as they rise in the column 


can be largely avoided by adding electrolytes or, preferably, surface-active agents. 


Résumé Des donneées sont presentees sur la dimension des bulles en colonnes a laide de quelques 
modéles. L’agglomeération des bulles fut photographice a des niveaux divers au-dessus des 
distributeurs de gaz. Les courbes de fréquence pour le diameétre idéal des bulles, a premiére 
approximation, sont symmetriques, a sommet unique et déterminées par le diamétre moyen, d,,,, 
et la déviation standard o. La dimension des pores filtrantes des distributuers n'a qu'une trés 
faible influence sur la dimension des bulles produites. Avec un debit gazeux croissant le diamétre 
moyen, d,,. des bulles s'accroit a un point fixe de la colonne jusqu’é ce qu'il atteint une valeur 
limite. Au fur que le niveau d’observation dans la colonne est élevé, toutes autres conditions 
restant égales, le diamétre moven des bulles, d,,. s'aceroit. Le diamétre de la colonne a une 
influence relativement faible sur la dimension des bulles. On peut éviter largement la coalescence 
des bulles ascendantes dans la colonne en ajoutant des électrolytes ou, de preference, des agents 


tensiouctifs. 
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KINLEITUNG 


BuaAseNsAULEN gehéren zu den dynamischen 


Zwei-Phasen-Systemen, bei denen cine Phase 
zerteilt in der anderen enthalten ist, wobei dieser 
Zustand durch 
beider Phasen zueinander aufrechterhalten wird, 
Bei den Blasensaulen steigt dabei cin Gasblasen- 


Kine 


gibt 


einen dauernden Gevenstrom 


strom in einer Fliissigkeitssaule auf. 


liber die mdglichen Systeme 


Sie sind fir die Verfahrenstechnik von 


Ubersicht 
Tabelle 1. 
Tabelle 1. Duynamische Zweiphasensysteme mit 


Phase 


einer zerteilten 





Phase 
serteilte 


suscmmenhdaugende Systembezcichnung 


Flussig fest Flissigkeitswirbelschicht 


gasformiyg fest Gaswirbelschicht 


fliissig fliissiy Tropfenstrom in 


Flissivkeiten 


fliissiy | gasformig | Blasensiiule 


gasfOrmig fliissig Tropfenregen im Gas 





grosser Bedeutung. Ejinmal lasst sich mit ihnen 
Raum 


Phasen 


auf kleinem eine grosse Grenztlache 


zwischen den erziclen, zum anderen 
bewirken die Partikel der dispergierten Phase 
durch ihre Bewegung cine Durchmischung des 


gesamten Mediums. 
Tabelle 1 
ihre 


Untersuchungen an den in der 


aufgefihrten konnen sich auf 


Erzeugung, ihre Dynamik und auf Austauschvor- 


Systemen 


ginge zwischen den beiden Phasen erstrecken. 


Es ist offensichtlich, dass dabei die geometrische 
Beschaffenheit der Elemente der verteilten Phase 
zentrale Stellung einnimmt. Sowohl thre 
Grosse als auch ihre Gestalt werden die Vorginge 


eine 
in den Systemen beeinflussen. Dieser Eintluss 
lasst sich in Wirbelschichten ohne Schwierigkeiten 
Die Partikel 


beliebiger Form und Grdéssenverteilung vorgeben 


untersuchen. festen konnen in 


werden, sie andern sich bei der Aufwirbelung 


nicht*. 


*Bei sehr kleinen Partikeln konnen Haftkrifte auftreten, 


sodass Teilchen in schwer kontrollierbarer Weise zu 
groésseren Elementen sich zusammenlagern, jedoch kann 


das als ein Sonderfall betrachtet werden. 


Ganz anders lieven die Verhaltnisse bei den 
Bei allen 
fluiden 


letzten drei Systemen der Tabelle 1, 
Kinrichtungen zur Dispergierung einer 
(strOmungsfahigen) Phase ist die absolute Grésse 
sowie dic Gréssenverteilung der entstehenden 
Teilchen abhangig vom Durchsatz. Dartiberhinaus 
aber kénnen die primar am Verteiler gebildeten 
Tropfen oder Blasen innerhalb des Tragermediums 
sich vereinigen oder zerfallen. 

Die vorliegende Arbeit soll am Beispiel von 
Luft-Wasser mit 
Abhangigkeit 


Blasengréssenverteilung vom Gasdurchsatz, von 


Blasensaulen des Systems 


porésen Gasverteilern die der 


der Héhe in der Saule und vom Siulendurchmesser 
aufzeigen. Dariberhinaus soll berichtet werden, 
Blasen durch 


Flissigkeit 


wie die Groéssenverteilung der 


Zusitze 


beeinflusst werden kann, wobei cine Homogenisie- 


obertlachenaktive zur 
rung der Blasengrésse in der Saule angestrebt 
wird. Es ist dabei nicht so sehr die Entstehung 
der Blasen 
vielmehr ihr Verhalten in der Saule Gegenstand 


unmittelbar am pordsen Boden als 
der Betrachtung, sodass die mitgeteilten Ergeb- 
nisse auch fiir Saulen mit anderen Gasverteilern 


gelten. 


LITERATURUBERSICHY 


Die Vorviinge in Blasenséiulen lassen sich zum 
Teil aus dem Aufstiegsmechnismus von einzelnen 
Blasen in Flissigkeiten erklaren. Die iiber solche 
Kinzelblasen erschiencne umfangreiche Literatur 
Siemes [1] und 
Hier soll daher auf eine solche 


wurde von zusammengestellt 


kritisch referiert. 
Zusammenstellung verzichtet werden. Aus einer 


Literaturtibersicht des gleichen Verfassers [2] 


uber die Entstehung von Gasblasen an Diisen und 


pordsen Platten seien nur die wichtigsten Arbeiten 


uber porése Verteiler zitiert. Gegenstand dieser 
Arbeiten Abhangigkeit der 
Blasengrésse von den Gas- und Flissigkeitseigen- 


schaften, von der Beschaffenheit der Gasvertciler 


war durchweg die 


und vom Durchsatz. 


Hateerstapr und Pratssxirz [8] hatten 
vefunden, dass die Blasengrésse bei sonst gleichen 
Bedingungen von der Art des verwendeten Gases 
In spateren Arbeiten anderer 


Als 


fir die Blasengrésse wichtige Fliissigkeitseigen- 


unabhangig ist. 
Autoren wurde dieses Ergebnis bestatigt. 





Ane. 1. Blasensiiule im Laminarbereich (Luftblasen in 


Wasser). 





Ann. 2. Blasensiiule im Turbulenzbereich (Luftbhasen 


in Wasser) 





Zur Blasenyrésse in Blasensiiulen 


schaft wurde von Hatsersrapr und Praussnrrz 
[3], Praussnrrz [4 5, 6] und Rupo.rn [7] die 
Oberflachenspannung, von ScuNURMANN [8, 9} 
Den 
Elektrolytzugaben zur Flissigkeit auf die Blasen- 
Kavursky und Turece [10] zur 
Erzeugung besonders feiner Blasen aus, 

Die Abhangigkeit 
Porositat der Gasverteiler wird von verschiedenen 
Forschern Nach 
ScHNURMANN [8], Motstap und Suvutman [11] 


die Viskositat herausgestelle. Einfluss von 


grésse nutzten 


der Blasengrésse von der 


unterschiedlich angegeben. 


sowie Hoveuron et al. [12] findet sich keine 
Abhangigkeit, nach Re vourn [7] und Pereruans 
[13] steigt die Blasengrésse mit der Porenweite 


an. Stemes und Borciers [14] zeigten, dass bei 
Luftblasen Wasser 
mittleren Porendurchmessers auf das Hundertfache 


in eine Vergrésserung des 


bei Metallfilterplatten nur eine Verdopplung des 


Blasendurclimessers zur Folge hat. Die Annahme 


einer Unabhangigkeit der Blasengrésse von der 


Porenweite seitens der erstgenannten Autoren 


lag somit nahe, da bei ihnen die untersuchten 


Porengréssenintervalle relativ klein waren. 
Mit 
Verteiler beobachteten verschiedene Autoren [7, 


wachsendem Gasdurchsatz durch den 
12, 13. 1 +} ein Anwachsen des Blasendurchmessers 
bei kleinen Durchsatzen bis zu einem Grenzwert, 
der sich bei weiterer Durchsatzsteigerung nicht 
mehr andert*. 

Das hydrodynamische Verhalter von Blasen 
siulen war Gegenstand Arbeiten von 
Suutman und Moustap [11], Verscnoor {15}, 
Stemes [1], Kéneen undStemes [16] und Borcuers 
{17}. Ks 
Zustande 
Gasdurchsatz. 


ciniger 


existieren danach zwei verschiedene 


solcher Saulen in Abhangigkeit vom 
Bei kleineren Durchsatzen steigen 
die Blasen ahnlich wie Einzelblasen auf, wobei 
die Flissigkeit relativ ruhig bleibt. Dieser Bereich 
wird als Steige- 
rung des Durchsatzes hat von einem recht genau 


Laminarbereich ” bezeichnet. 
erkennbaren Wert an cine starke Beunruhigung 
zur Yolge, die wieder auf die 


Blasen Ks 


sich dauernd in Zusammensetzung und Bewegung 


der Flissigkeit 


aufsteigenden zuruckwirkt. bilden 


wechselnde “ Turbulenzballen “ aus, die das 


*Das gilt jedoch nur innerhalb des unten definierten 
Laminarbereiches, 


79 


Innere der Saule stark verwirbeln, so dass ein 


Gas-Flissigkeits-Gemisch mit stets wechselnder 
Dieser Durchsatzbereich 
Die Abb. 1 


veranschaulichen diese beiden Stro6mungszustandei 


(1) geklart 


Phasengrenze entsteht. 


ist der “* Turbulenzbereich.”’ und 2 


deren Entstehung von Stemes 
wurde. 

Die 
Blasensaulen ist am leichtesten durch die Aus- 
Methode 
HALBERSTADT und 
Die an Filterplatten 
entstehenden Blasen sind nicht gleich gross, ihre 


exakte Erfassung der Blasengrésse in 


wertung von Fotos mdglich. Diese 


wurde wohl zuerst§ von 


Praussnirz [3] benutzt. 


Durchmesser tiberde*ken ein ganzes Spektrum. 
Siemes [1] zeigte, dass sich dieses Spektrum 
niherungsweise durch eine Normalverteilung 
wiedergeben lasst*. Eine erste cingehende Unter- 
suchung der Blasengréssenverteilung in einer 
Saule von 92 mm. lichter Weite mit porésem Boden 
verOffentlichten Siemes und Borcners [14], die 
Methoden 
Da 


Auswertung von Fotos aus dem Turbulenzbereich 


sich bei der Auswertung auf von 


Siemes und GUnruer [18] stitzten. eime 
kaum mdglich ist. weil hier die Blasen sehr stark 
und unregelmissig deformiert sind und sich die 
Teil 


») 


Blasengrenzen = zum nicht  identifizieren 


Abb. 


Mitteilung auf den Laminarbereich. 


diese 
Die Vertei- 
gekennzeichnet 


lassen (sieche beschrankt = sich 


lungsfunktionen werden durch 
Mittelwert 


Diese Gréssen werden fiir Metalltilterplatten mit 


den und die Standardabweichung. 
verschiedenen Porenweiten in Abhangigkeit vom 
Als Flussigkeiten 


wassrige 


Gasdurchsatz wiedergegeben. 


wurden Wasser sowie Losungen von 


Kaliumathylxanthat und Aluminiumsulfat 
verwendet. 

et al, 
Unkenntnis friherer Arbeiten, Untersuchungen 
Saule 


Boden durch 


HowuGuron [12] fuhrten, offenbar in 
78 mm 
Filter- 
platten nur teilweise begast wurde. Sie teilen die 


der Blasengrésse in einer von 


Durchmesser durch, deren 
fir verschiedene 


mittleren Blasendurchmesser 


Fliissigkeiten und Filtermaterialien in Abhangig- 

*Bei an einzelnen Diisen entstehenden Blasenschwirmen 
stimmt das und 
GUntruer [18]); hier kann es zweckmiassig sein, statt 
der Hiaufigkeitsverteilung fiir die Durchmesser die fiir 
Blasenvolumen oder die Blasenobertlachen zu betrachten. 


im allgemeinen nicht mehr (Stemes 





W. Stemes und EK. Borcurrs 


keit vom Durchsatz mit, wobei sie sich ebenfalls 
auf den Laminarbereich beschrianken. 


VERSUCHSDURCHFUHRUNG | 
AUSWERTUNG 


a. Apparatur 


Abb. 38 


Apparatur wieder. 


gibt die wichtigsten Elemente der 


Die Blasensiulen wurden 


(1) 


2m hohen Aecrylglasrohren erzeugt. Die 
































’ 


») 


3. (1) 2 


Manometer, (3) Gasulr, (4) Regulierventil. 


Ans. Versuchapparatur : \crylglasrohr, ( 


kilter- 
platten uber den gesamten Bodenquerschnitt der 


Gasverteilung erfolgte durch keramiusc! 


Saulen. Uber den Eintluss der Filte rporengrosse 


auf den Blasendurchmesser wurde bereits friiher 
[14] berichtet. Hier sollen nur Versuche herange- 
Platten 
mittleren 


zogen werden, bei denen von 20mm 


Dicke 
von 25 bis 35 verwendet wurden. 
Platte 
zufihrung und Flussigkeitsablass. 


mit emem Porendurchmesser 
Unter der 
mit Laft- 
Druckdif- 


Gaskammer 
Di 


Aussenraum 


befand sich eine 


ferenz zwischen dem der 


Druckluftzufihrung am unteren Saulenende 
wurde mittels eines Quecksilbermanometers (2) 
gemessen. Die durchyesetzten 
Pressluft erfolgte (3). 


Davor war noch ein Regulierventil (4) geschaltet. 


Messung der 


mittels einer Gasulr 


Der Gasdurchsatz wird im folgenden stets auf 


die Einheit des Saulenquerschnitts bezogen 


80 


angegeben, damit die Ergebnisse fiir verschieden 


unmittelbar vergleichbar — sind. 
Zahlenwerte 


umygerechnet. 


weite Saulen 


Die 
Atmospharendruck 


sind auf 


Bei 


mit 


angegebenen 
den 
Messungen) wurden  Acrylglasrohre cinem 
inneren Durchmesser von 42, 92, 192 und 292 mm 
Die Siiule 
betrug ohne Besagung jeweils 90 em. 

Neben Wasser 


wissrige Losungen untersucht, und zwar Wasser 
] 


verwendet. Flussigkeitshéhe in der 


destilliertem wurden zwei 


mit Elektrolytzusatz | Aluminumesulfat 
Wasser mit 


Substanz (0.5 ¢ Kaliumathyixanthat pro Liter). 


uw 


pro 
Liter) und einer obertlichenaktiven 
Der Wahl dieser Lésungen lag die Beobachtung 
Zusiitze der 
Blasen- 


\uswirkungen 


zugrunde, class schon voringc 


genannten Substanzen Anderungen des 


verhaltens zur Folge haben, deren 
auf das Verhalten der 
Modellfallen untersucht werden sollte. 

Dic mit Leica auf Leitz 
Wechselschlitten (Naheimstellgerat ) 
(f $. 5) 

Objektiv etwa 40 45 em von den Saulen entfernt 
Die Kamera Elektronen 
blitzverat svnchronisiert, das sich in Kamerarich- 
Licht 
durch mehrere Bogen Transparentpapier diffus 
gemacht Diese 
konnte in der Hohe verstellt werden. 
Fehler durch die Zufalligkeit des Bildausschnittes 
jeder Messung zwei 


Blasensaulen an dliesen 


kotos wurden cmer 


mit Focotar 


jeml: aufgenommen. wobei das 


war. war mit emem 


tung hinter der Saule befand und dessen 


wurde. Aufnahmeeinrichtung 
Um vrobe 
wurden bei 


zuesvermeiden, 


Fotos in etwa 20sec Abstand aufgenommen. 


bh. 


Auswertung der Messergebnisse 


Dic 


unter Verwendung eines 


erfolyte 
Dic 
projizierten Negative zeigten die Blasen in etwa 
ie 

der 


Auswertung des Filmmaterials 


Mikrotilmygerates. 
t-facher natirlicher Grosse. Die Bestimmung 


realen Abmessungen erfolgte tuber eine 


mitfotogratierte Vergleichsskala. \usgemessen 
wurde fiir jede Blase der maximale Durchmesser. 
Die Blasen wurden nach diesen Durchmessern in 
Da die 


fur 


verschiedene Groéssenklassen cingeordnet. 


Hautigkeitsverteilungen aber nicht clie 


maximalen, sondern fiir die ideellen* Durchmesser 


*Der ideelle Durchinesser einer Blase ist der Durchmesser 


der ihr volumengleichen Kugel. 





Zur Blasengrésse 


Annahme 
Nach SiemMes 
{1} und Stemes und GOwruer |18} wurden die 


sollten, 


lit 


liber die Gestalt der Blasen nétig. 


werden war eine 


Blasen als Kugeln (Durchmesser unter 1 mm) und 
Die 


gemessenen grossen Durchmessern der Ellipsoid- 


Rotationsellipsoide —idealisiert. m—Ssodden 
blasen gehérenden kleinen Durchmesser wurden 
als eindeutige Funktion des Maximaldurchmessers 
und der Fliissigkeit betrachet. In Vergleichsmes- 
Funktionen 


50 Blasen gleichen Maximal- 


sungen wurden diese bestimmt, 


indem fir je etwa 
durchmessers die kleinen Durchmesser ausgemes- 
sen und gemittelt wurden. Aus den gemessenen 


Maximaldurchmessern und den thnen zugeord- 


neten kleinen Durchmessern liess sich das Blasen- 
ideclle 


Kinzelheiten 


damit der Durchmesser 
Wevgen 
Auswertemethode sei auf Stemes und GUNTHER 

7 Sie 
Form auch Hor GHTON 
In Tabelle 2 sind fiir Wasser 


die verwendeten Groéssenklassen fur die Maximal- 


volumen und 


bestimmen., der clieser 


[18] sowie auf Borcuers [17] verwiesen. 


wurde in éhnilicher von 


et al. |12| verwendet. 


durchmesser und die ihnen zugeordneten Werte 


der kleinen und ideellen Durchmesser wieder- 


geveben. Fir die wiassrigen Loésungen sind die 


entsprechenden Werte davon nicht sehr verschie- 


Tabelle 2. 


Durchmesser von Lufthlasen in Wasser 





ike ‘ lhe ’ 


Durchmesser 


mittleres 


grosster hleines 


Durchmesses Durchmesses d 


(trim) (rrimn) (rim) 


~ we w tt 





Die auf den ideellen Durchmesser d bezogenen 


Hautigkeitswerte A, fur die einzelnen Gréssen- 


klassen / mit den mittleren Klassendurchmessern 
Klassenbreiten Ad; 
Blasenzahlen N, 


d;, den und den ihnen 


zukommenden ergeben sich zu 


in Blasensiulen 


(1) 


Sie stellen in guter Naherung Gauss’sche Vertei- 
lungen dar, wie die Abb. 4 und 5 an Beispielen 


—— Sdulenhone = 30°5 cm 
——— Sdulenhone = 80*5cm 








Blosenhoufighet, 

















Ann. 4 fiir Luftblasen in 


Wasser bei Ve 


Haufigkeitsverteilung 
0.38 ems in einer Siule von 92 mm 
Weite (Aufnahmehohe : 
platte). 


300m oberhalb der Boden- 


durch 
den mittleren Durchmesser d,, und die Standard- 


zeigen. Solehe Verteilungen lassen sich 


abweichung o kennzeichnen. d,, kann in guter 
Naherung nach 

ZN, d; 

ZN; 
berechnet werden. Fur das Quadrat von 
dann : 

2 N; (d,, 
ZN 


d; )? 


: 
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Luftblasen 
J ems in einer Sdule von 42 mm 

SOem oberhalb der Boden 
platte) 


Ann. 5 Hliutivkeitsverteilung = fiir in 


” 


Wasser bei | 
Weite 


W 
( Nufiiaborne doh 


VERSUCHSERGEBNISSE 
Die 
Durchsatz 


a. Abhingigkeit der Blasengrisse vom 


Den Eintluss des Durchsatzes auf die Blasen- 
grésse zeigt die Abb. 6. 
d,. fir die 
Abhangigkeit 
Die Blasendurchmesser ergaben sich aus 
wo 


Hier sind die mittleren 


Durchmesser untersuchten Flissig- 
keiten 
tragen. 


Aufnahmen, 


in vom Durchsatz aufge- 


die in emer Zone von - Sem 


oberhalb der Bodenplatte in der 92 mm-Saule 
In 


bereichen andert sich dic Blasengrésse schr stark, 


gemacht wurden. den unteren Durchsatz 


Der mittlere Blasendurchmesser verdoppelt sich 
in einem Durchsatzintervall von 0,01 — 2.0 cm sec, 
Die grésste Zunahme erfolgt bei Wasser, namlich 
von Imm auf etwa 25mm, die geringste 
Vergrésserung der d,-Werte beobachtet man bei 
den obertlichenaktiven Lésungen, wahrend die 


Elektroly thisungen cme Mittelstellung cinnehmen. 


sZ 











%. cm sec 

Mittlerer Blasendurchmesser als Funktion 
(Laft 
Lasungen) fiir eine Siule von 02 num Weite (Aufnahme- 
Socom oberhalb der Bodenplatte) 


Ann. 6. 


des Durechsatzes in Wasser und wiissrigen 


Beim weiteren Anwachsen des Durchsatzes wird 
d,, bleibt 
oberhalb cines Durchsatzes von 3 em/see nahezu 
konstant. Solche 


fiir 


nur noch unwesentlich grésser und 


Abhingigkeiten ergeben sich 


auch andere Filterplatten, andere Saulen, 
Sic 


mehrfach 


sowie fir verschiedene Hihen in der Siule. 


Abnlicher Form bereits 


publiziert [12, 14, 17), 


wurden in 


+ 














om sec” 


K 


Ann’ 7. 
Durchsatz fiir die in Abb. 6 angegebenen Bedingungen, 


Abhidingigkeit der Standardabweichung vom 


Abb. 7 
Abweichungen fir die 92 mm-Saule in 30 + 5 em 


kleinen 


bringt die mittleren quadratischen 


Die Werte steigen wie die d,,-Werte bei 


Durchsaitzen bis zu 2em/sec stark an, 


um dann fast konstant zu bleiben oder doch nur 





Zur Blasengrésse in Blasensiiulen 


Die 
Werte fiir Wasser liegen wieder tiber denen der 
Auch der 
Verlauf der Kurven 


noch einen geringen Anstieg aufzuweisen. 


hier wird qualitative 
Filter- 


Beobachtung- 


Lésungen. 


bei Variation von 
platten, Saulendurchmesser und 
Auch bei der Auftragung 
findet 


dieser Anstieg der Streuungsbreite bei kleinen 


shéhe nicht geindert. 


der relativen Abweichungen  o/d,, sich 
Durchsatzen. 
Die Deutung 
Betrachtung des Bildungsmechanismus der Blasen 
Filterplatten Bei 


Durchsatzen wenige 


der Ergebnisse kann an die 


an den anschliessen. sehr 
kleinen 
aktiv. 


flussen 


sind nur Poren 


Die kleinen aufsteigenden Blasen beein- 
sich noch nicht, eine Vereinigung zu 
grésseren findet kaum statt, daher sind sowohl 
die d,, Werte klein. Mit steigen- 


dem Durchsatz wird der Abstand der gasenden 


als auch die o 


Poren kleiner, bis er der Grésse der entstehenden 


Blasen entspricht. Die Vereinigung kleiner Blasen 


zu grésseren vollzieht sich immer starker, wobei 
\ufstiegsgeschwindigkeit = der 


Blasen 


schliesslich durch den Aufstieg aus der Bodenzone 


gleichzeitig die 


resultierended anwichst, und — sie 


einer weiteren Vergrésserung entzogen werden. 


Ann. &. 


fiir verschieden 


Y, 


Mittlerer Blasendurchmesser als Funktion des Durchsatzes (Luft 
weite Siulen (Aufnahmehdhe : 


Art 
und 


Es stellt sich dann offensichtlich eine 
Gleichgewichtszustand Bildung 


Aufstieg von Blasen ein, der zu den mittleren 


zwischen 


Blasendurchmessern in den grésseren Durchsatz- 
Die d Werte sind daher auch 
Bereich nahezu konstant. Abhnliches 
Werte. Hat Abstand 


einzelnen gasenden Poren Blasengrésse erreicht, 


bereichen filhrt. 


m 
in diesem 
gilt fiir die o der der 
wird das Blasenspektrum durch teilweise Vers- 
chmelzung kleiner Blaschen breiter, was zu einer 
Werte fiihrt. Ist der 


Bildung und 


starken Zunahme der o 


Gleichgewichtszustand zwischen 
Aufstieg 


wesentlich. 


erreicht, andert sich o nicht mehr 
Auf die Unterschiede bei den unter- 
suchten Fluissigkeiten wird weiter unten eingegan- 


gen. 


hb. Die Abhingigkeit des Blasendurchmessers vom 
Sdulendurchmesser 

Die oben fiir die 92 mm — Saule_ beschriebene 
Abhangivkeit 
beobachtet man in Ahnlicher Form auch bei den 


der Blasengrésse vom Durchsatz 


anderen Saulen. Die dabei auftretenden d,, — und 
o Werte zeigen jedoch einen geringen Einfluss 
Abb. 8 zeigt das fiir 


des Saulendurchmessers. 








cm sec 


in Wasser) 
s0em oberhalb der Bodenplatte). 
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Wasser bei ciner Messhéhe von 30 + 5¢m. zu den grésseren stirkere Streuungen der Werte. 


Hier sind die d,— Werte in Abhangigkeit vom Die Standardabweichungen — liegen, abyveschen 


Durchsatz fiir verschieden weite Saulen aufge- vom Bereich kleiner Durchsatze (1 2 om sec), 
tragen. Abb. 9 stellt die entsprechend no - Werte fur alle Saulen gut beieimander. 
dar. Bei den Kaliumiathylxanthat-Lésungen findet 
sich der Einthuss des Saulendurchmessers auf die 
42mm Sdule d,, Werte nicht mehr so ausgepragt wie bei 
ne Wasser, wie aus der Abb. 10 hervogeht. Hier 
© 292 0m Séule fallen nur die d,— Werte fiir die 42 mm — Saule 
aus der Kurvenschar heraus. Gleiches gilt fiir die 
Standardabweichung (Abb. 11). 
Die Elektrolytlisungen nehmen auch hier 


wieder eine Mittelstellune ein. Der Einfluss 


Ann. 0. Wie Abb. 7, jedoch fiir die in Abb. & anges 


benen Bedingunyen 


Die mittleren Blasendurchmesser nehmen unter 
sonst gleichen Bedingungen mit wachsendem 
Siulendurchmesser ab, sofern der Saulendurch 
messer unter 200mm bleibt. Fir Saulen mit 
einem Durchmesser tiber 200mm _ bleibt 4, 
unabhangig vom Saulendurchmesser. Die sehr +, jedoch fiir die 


enge Siule (42 mm 4) zeigt ferner im Gevensatz u nen Bedingungen 


cm sec 


Ane. 10. Wie Abb. 8. doch fiir Luft in wissriger Kaliumiithvixanthat Lisung. 
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auf Randeffekte 
zurikzufiihren. Je Saule ist. 
starker machen sich die Randwirkungen bemerk- 


des Saéulendurchmessers — ist 


enger die desto 
schwankenden 
Saule 


Verhaltnis 


den erheblich 
der 

Niederschlag findet. Errst 
der Querschnitte von Randzone und stérungsfreier 
Innenzone einen bestimmten Wert unterschritten 
hat, 
vernachlassigbar. 
192 mm — Saule 
Kurvenverlauf nahezu identisch mit 
Saule. Als Richtwert kann 
dass ein Randeintluss der Saule nicht mehr zu 
beriicksichtigen ist. wenn der Saéulendurchmesser 


bar, was auch in 


d Werten 


m 


bei 42 mm seinen 


wenn das 


Einfluss des Siulendurchmessers 
Dieser Wert der 
erreicht) zu 


wird der 


scheint bei 
Hier ist 


dem 


sein. der 
der 


202 mm gelten, 


50 mal grisser ist als der mittlere Durchmesser 
der Blaser ° 


«. Die Abhidngigkeit des Blasendurchmessers von 
der Hohe 


schafti 7 


in der Sdule und den Fliissigkeitseigen- 


Den Eintluss der Hohe in der Saule (das ist die 
Hihe. in der die Saule fotografiert wurd) auf die 
d,,- Werte geben die Abb. 12. 13, 14, 
fir Wasser und verschiedene Siulendurchmesser 


und die Abb. 16 und 17 fiir die Lésungen in der 


und 15 


02 mm — Siule wieder. 


Sdulenhohe 30+ Scm 
Saulenndhe 60+ Scm 
Sdulenhohe BO* Scm 





3 
4  ocmsec 


Ase. 12. Mittlerer Blasendurchmesser als Funktion 


Durchsatzes (Luft in Wasser) fiir verschiedene 


des 


\ufnahmehohen (Siulenweite : 42 mm). 


Durchsitzen ist = die 


Blasendurch- 


Bei sehr kleinen 
Hoéhenabhangigkeit des mittleren 


messer so gering, dass man sic naherungsweise 
Das liegt daran, dass die 
entfernt 


vernachlassigen kann. 


Blasen so weit voneinander 


bleiben, 








Saulenhohe 30°5cm 
Saulenhohe 6015cm 
© Sdulenhdhe 802 5cm 


cm sec 


Wie Abb. 12, doch Siulenweite 92 mm. 


Saulenhohe 30° 5cm 
Séulenhohe 60* 5 cm 
Sdule hohe BO* Scm 


Ann. 14. Wie Abb. 12, doch Séiulenweite 192 mm 


Soulenhuh 
Soulenh 


Wie Abb. 12. doch Saiulenweite 292 mm. 


¢ Soulenhohe 80+ Scm 
Saulennéhe 60+ S5cm 
Sdulenhdhe 30* Scm 











cm sec 
Wie Abb. 14, doch Luft in’ wiissriger 
A\luminiumsulfat-Losung. 


Ann. 16. 
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die Elektrolytlésungen. Hier kann man wahr- 


Sduienhohe BO Sem scheinlich eine Ausbildung — elektrostatischer 
© Sédulenhohe 60° Serr 
« Sdulenhohe 30° Sem 

a eingehend von Corux und Neumann [19] und 


SCHNURMANN [20] diskutiert worden ist. Infolge 


Doppelschichten um die Blasen annehmen, was 








der elektrostatischen Aufladung stossen sich die 
Blasen ab, so dass eine Vereinigung zu grésseren 





2 3 Blasen erschwert wird. 


& ome Mit) wachsendem  Siulendurchmesser — wird 


bei Wasser der Héheneinfluss auf den Blasen- 
Wie Abb. 14. doch Luft in wiissriger 


syd durchmesser kleiner. Vermutlich tindet im beson- 
Kaliumaethylxanthat-Losung. 


ders stark verwirbelten Grenzbereich der inneren 
Aufwartsstromung und der wandnahen Abwirts- 
. ; . . str uo ‘ime bevorz ote ase -ereinie wv 
dass sie nicht zusammenlaufen kénnen, Erst StOmung ein vorzugte Bla caverns 
ems dew Duschasta ebwe 09 ~@.4 caniece crveicht statt. dieser Bereich ist aber bei engen Saulen 
hat. treten bei Wasser erkennbare Differenzen 


zwischen den d,,— Werten der verschiedenen 


relativ am gréssten. Wahrend bei Wasser in der 

92 mm —- Siule der mittlere Durchmesser fiir eine 
A ‘s ¥ > > . g . os oy Bet > vn 

Aufnahmehdhen auf*. Diese Vergrésserung des Aufstiegsstreke von 50emumea. 25 per cent, 


mittieven Blasendurchencesers mit der Hithe anwichst, was einer Verdopplung des mittleren 


Blasenvolumens entspricht, betragt die Zunahme 
der d Werte fiir die gleiche Strecke in der 


geht auf eine Vereinigung der Blasen beim 
Aufstieg zuriick, sie ist stark abhangig von der 
Art der Flissigkeit. Man beobachtet bei den 
oberflachenaktiven  Lésungen praktiseh keine Anerkennung Dem Leiter des Institutes tir Technische 
Verinderung des Blasendurchmessers (Abb. 17), Chemie, Herrn Professor Dr. phil, H. Kounet danken wir 
fiir sein wohlwollendes Interesse und die stete Forderung 
der Arbeit. Wir danken auch der Max-Buchner-For- 
; schungsstiftung fiir die Bereitstellung von Mitteln zur 
gen wieder cine Mittelstellung ein. Beachaffung der Apparaturen. 

‘ Die erhohte§f Homogenitat der Blasen in den 

oberflachenaktiven Losungen lasst sich mit der SYMBOLVERZEICHNIS 

Bildung von stabilisierend wirkenden Blasen- 


292 mm — Saiule nur etwa 10-15 per cent. 


was auch fiir die o Werte gilt. Wie aus der 
Abb. 16 hervorgeht, nehmen die Elektrolytlésun- 


* ideeller Durchmesser 
grenzschichten erklaren. Ahnliches gilt auch fur mittlerer Klassendurchmesser 
mittlerer Blasendurchmesser 


*Die Vergrésserung von d,, auf Grund des nach oben Hiutigkeitswerte fir die cinzelnen Grossenklassen / 


abnehmenden hydrostatischen Druckes betrigt etwa Anzahl der ausgezihiten Blasen in der klasse / 


1.5%, fallt daher hier in den Fehlerbereich der Bestim- : Gasdurchsatz pro Zeit) und Flicheneinheit 
mungsmethode. ‘ Standardabweichung 
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On optimum cross-current extraction 


R. Ari. D. F. Repp and N. R. Amunpson 


Department of Chemical Engineering, University of Minnesota 
(Received 1 October 1959) 


Abstract—The notion of dynamic programming is here applied to the problem of allocating 
an extracting solvent to the various stages of a cross-current extraction unit. The basic ideas of 
optimizing a multi-stage process by the techniques of dynamic programming are first explained 
and then shown in detail in an example of extraction by an immiscible solvent. This problem 
is fully solved either when a value is ascribed to the solvent or when it is available in a fixed 


total quantity, Finally it is shown how the case of miscible solvents can be reduced to a similar 


computation, 


Résumé—Les auteurs appliquent la notion de programmation dynamique au probléme de 
lalimentation du solvant d’extraction aux divers étages dune unite dextraction a contre- 
courant. Ls expliquent d'abord comment rendre optimum un procédé a étages multiples par les 
techniques de programmation dynamique et montrent en detail un exemple dextraction par 
solvant non miscible. Le problonne est entiérement résolu pour une fraction préevue de solvant 


comme pour la quantité totale fixée. Le cas des solvants miscibles peut étre ramené a un calcul 


analogue. 


Zusammenfassung Das Verfahren der dvnamischen Programmiecrung ist hier auf das Problem 
der Zuteilung des Extraktionsmittels auf die einzelnen Stufen einer Kreuzsrtomextraktions- 
apparatur angewandt. Die Grundsitze der Programmierung eines Vietstufenprozesses mit Hilfe 
der dynamischen Programmierung werden erliutert und dann im einzelnen an cinem Beispiel 
der Extraktion durch eine unmischbare Flissigkeit gezeigt. Das Problem ist vollstandig gelist 
fiir den Fall, dass man dem Loésungsmittel einen Wert zuschreibt oder dass es in einer bestimmten 
Gesamtmenge verfiigbar ist. Albschliessend wird noch gezeigt, wie man den Fall des mischbaren 


Losungsmittels auf eine cinfache Rechnung zurtickfihren kann, 


Tur Norion or Dynamec guided by a gradient technique. What is needed 
PROGRAMMING is an algorithm which reduces the dimensionality 
Tur problem of choosing the conditions in a of the problem by replacing it by a sequence of 


multi-stage process to achieve an operation in problems with fewer variables. Such an algorithm 


some sense optimum has always been of the first 1s provided by the notions of dynamic programm. 


importance in the chemical industry. In earlier ing which have been put forward by Brian 


vears it was part of the craft of the plant operator, and his co-workers. This work, begun some ten 


but with more widespread use of automatic Years ago, and first published in a Rand Corpora- 


control there has been an increasing demand for tion monograph in 1953 has received brilliant 


methods of design and calculation that will give exposition in numerous papers and in the book 


optimum conditions. For processes with many [1]. The elementary example of multi-stage, 
operating variables the standard method of  cross-current extraction which is given here will, 
setting the derivatives of the profit with respect 
equal to difficulties of the method, but the underlying 


it is hoped, bring out some of the virtues and 


of all the variables simultaneously 


zero rapidly leads to equations of intolerable Notions are sulliciently beautiful to merit a more 


complexity. Direct numerical calculations in abstract presentation. 


search of the maximum are laborious even when Let the process consist of N stages and the 
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process stream be characterized by a vector p. 
P, denotes the state of the feed stream to the 
first stage, p, the state of the product of the first 
stage, which is also that of the feed to the second, 
and py is the state of the final product stream. 
Each stage transforms the state of the stream in 
a WAY which depends on the values of the operat- 
If these operating 
Gy... -y. 
where the elements of q, are the values of the 


th 
i 


ing variables of the stage. 


variables are denoted by vectors 


operating variables for the stage, we may 
P, = T(Pnas Ga) 
P,, 7... t the 
the of the stream. 
related to py by the combined transformation 
Py Ty Ty ,..-Te7, Be If M 
elements in each q, there are thus a total of MN 
to be 


profit 


write, or more briefly, 


transformation 


Thus py, 


denote 


state process is 


there are 


variables whose values have 
Our to 


function P (py) of the state of the output stream. 


operating 
chosen. aim is maximize some 
As it stands the problem is thus to assign values 
to the MN operating variables and with even 
moderately large M and N this is a formidable 


The 


which 


choice 
the 
will be called the optimal N-stage 


business, correct of the operating 


variables achieves maximum 2 (py) 
policy with 
respect to the feed state py. for when this choice 


of q has been made py, Ty Ty ,...7, De isa 
function of the feed state py only. We may write 


Max P (py) 
4, 


ly (Po) (1) 


The N-stage process may now be thought of as 
a single first stage followed by an (NV 1 )-stave 
In the first stage the feed stream state 
T, Py and this 
the 


operating 


process, 
is transformed from py to p, 
the feed to the (N process ; 
transformation 7, depends on the 
variables q,. Certainly if the whole process is 


IS 1 )-stage 


operating with the optimal N-stage policy with 
respect to the feed state py. the last (NV 1) 
the optimal 
(N 1)-stage policy with respect to their feed 
P;- Any the 


variables q, than the values appropriate to the 


staves must be operating with 


other assignment of operating 


N-stage optimal policy will decrease the profit 


function, even though the last N 1 stages are 


of 


89 


(N 1)-stave with 


Thus 


using the optimal policy 


respect to their feed. 


q,)| 


fx (Po) Max fy 1 [7 (Po; 
a, 


This equation can be used to get Ts when fy , IS 
known, and moreover 


Max P |T (py; 


4 i 


/, (Po) q,) | (3) 
In cach case the maximization only involves the 
M variables of q, and we have replaced a problem 
in the simultaneous variation of MN variables 
by a sequence of N problems in the simultaneous 
variation of M variables. a much more tractable 
situation. 

of the method be 
Firstly, in calculating the sequence of 
functions /,; (p). fo (Pp) 


tion that has been accumulated is useful. viving 


Several advantages mia 
noticed, 


fy (p) all the informa- 


optimal policies for all numbers of stages and for 
Secondly, if. 
in many physical processes, the variables q are 
bounded the 
T (p; q) are continuous in p and q and P (p) is 


continuous in p then the maximum will certainly 


various feed states. as is the case 


above and below, transformations 


exist. Thirdly, the method is ideally suited to 


modern high speed computers, for they work 


most happily with a standard, oft-repeated 
algorithm such as is enshrined in equation (2), 
Moreover they use a direct search technique which 
is made easier by the restrictions on physical 
variables, in contrast with the analytical techni- 
ques which do not take kindly to restrictions in 
the of Ditliculties if 
the fy (p,) has a number of local maxima for 
q; it difficult to 
devise a search routine that will ensure that the 
of 


exhaust the storage space (even in the largest of 


form inequalities. arise 


admissible for is sometimes 


vreatest these is found. It is also easy to 


computers) if p has many components, 
In 


ponents of the vector p would be the temperature, 


a chemical engineering context the com- 


concentrations, flow rate, etc. required to pre- 
the of the The 
operating variables q might be the size of the stage, 


seribe state process stream. 
its operating temperature or temperatures, the 
amount of solvent used or flocculent added and 


so on, 
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The method outlined can be easily modified 
to take into account the cost of operation of each 


stage. Let the return from the process be the 


profit function P (py) diminished by the cost of 
N 


operation 2 C (q,) of each stage. Again an 
1 


optimum N-stage policy with respect to py will 
vield a maximum, 


Max [P (py) 


Sy (Po) 


\ 
XC (q,)} (4) 
1 


and an argument entirely analogous to the 


preceding gives 


Max { fy, [1 (Pos ay)] — C(ay)} (3) 


ty (Po) 
with 


A, (Bo) = Max { P [7 (Pos a) ] 


C(q,)} (6) 


In both the functional 2) 
and (5) have been obtained by the application 
This 


has been stated by Bettman [1] in the following 


cases equations 


of the so-called principle of optimality. 
terms : 


an optimal policy has the property that 


whatever the initial state and initial decision 
may be, the remaining decisions must constitute 
an optimal policy with regard to the state result- 


ing from first decision. 


N-STAGE 
EXTRACTION 


CROSS-CURKRENT 


We will consider the simplest of extraction 
The 


solvent passes through the sequence of stages at 


processes with immiscible solvents. first 
a flow rate q- The state of this solvent is specitied 
by a single variable x the concentration of solute 
(x plays the role of p, a vector with only one 
the of the 


paragraph). Through each stage there passes a 


element, in presentation previous 
stream of extracting solvent, which for convenience 
will be referred to as the wash water, with flow 
rate w,, to the n™ stage (w plays the role of q in 
the previous discussion). The wash water leaves 
the stage with solute concentration y,. which 
will be taken to be in equilibrium with the con- 


centration z, in the first solvent. The total wash 
\ 

2'w,, and the total amount 
i \ 

of solute extracted is Q, 


water used is Wy 


q (\%q ~ yy). 


a" ome 
ae, Vy, 


Rupp 


4 


« 


and N. 


0 


R. AmunNpson 
Two related optimum problems may now be 
stated. 


1. With a fixed total resource Wy of wash 
water we want to distribute it between the 
stages to maximize the total amount ex- 
tracted Qy. 

If the cost of providing and pumping the 
wash water is a fraction A of the value of 
the extract we wish to assign the w,, so that 


the profit Q\y — A Wy is maximum, 


Clearly in the first problem it is necessary to 
put an upper limit on the total amount of wash 
water else it would be possible to extract all the 
solute by using an infinite amount of wash water. 


©) @ 
Ose! ths, 


Yi 2 Yn 


Schematic diagram of cocurrent extraction. 


ss 
Rn-' 


Fic. 1. 


The second problem is more in line with the 
economics of the situation and here there is no 
need of a bound on Wy for the profit would 
become negative if Wy were too large. We shall 
show that the solution of the second problem 
vields also the solution of the first and is in fact 
the easiest method of solving it. Readers familiar 
with constrained maxima will recognize A as the 
Lagrange multiplier and will find an exposition 
of its use in the following references [1-3]. 
In the first problem suppose that ry = @ and 
Wy = 6 then when the optimum policy has been 
found the maximum extract Qy is a function of 
both a and b, 
(7) 


fy (a,b) = Max Q, 


ey) 


\ 
J 


b. (We note that @, 


p> mw, 


) 
increases monotonically 


when #, = @ and 


with Wy and so the 


‘: b). The 


optimal policy will not give Dw, - 


principle of optimality now gives 
Sy (a,b) — Max! fy , (4b — w,) + gla —a,)} (8) 


for if w, has been used in the first stage there is 
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left to be distributed 


1) stages. 


only an amount b — w, 
among the remaining (N This is a 
straight forward enough problem but it involves 
calculating a function of two variables, and this, 
in terms of machine storage, is expensive. 

In the second problem the maximum of the 
profit function Q, AW, is a function of the 
inlet concentration ry only when A has been fixed, 


We may write 


fy,, (a) = Max (Qy — A Wy) (9) 


ys 
and the principle of optimality gives 

axif ' 
Max} fy—3.a(4y) 


u 


ty (a) q(a —a,) — Aw, } (10) 
Here a function of only one variable fy , (a) has 
to be computed and stored during the sequence 
of calculations with given A. When the optimum 
policy has been found Wy. the total wash water, 
can be calculated and this will depend on the 


Wy Wy (A). When 


such caleulations have been done this function 


value of A, Le. several 
has been calculated and the value of A for which 
Wy (A) 
for this value of A gives the solution of the first 


b may be found. The optimal policy 


problem. In terms of machine storage and time 
this sequence of calculations is much cheaper 
than the direct solution of the previous problem. 
Two slight generalizations may be incorporated 
in these problems. In the first place we might 
wish to include the cost of providing and main- 
taining each stage; the profit function would 


then be @, AW, 


as a fraction of the value of the extract. 


uN where pu is this new 


cost 


Secondly, it might be that a different use should 


be made of, and hence a different value ascribed 
to. the extracts from each stage according to 
We might 


thus put a premium on obtaining higher concen- 


the concentration of solute emerging. 


trations in the extract and the value of the extract 


uM ”" Yn q { vy, 1 


might be zero for y,, 


r,) would be R (y,,) w,, y,. BR (y,) 
y* say, meaning that the 
extract is valueless if produced at lower concentra- 
These 
porated in the model to be described and further 
the final 
the 


tion than y*. modifications will be incor- 


generalizations will be discussed in 


paragraph, but we will first set up basic 


equations. 


3. THe EQuarions 


A mass balance on the solute around the n"' 


unit gives 


Wy Yn q (%-1 iy) 


and since q ts constant it is convenient to define a 


dimensionless tlow rate 


Then 


7 (11) 


Un Yn Un 1 


Since equilibrium is assumed in each stage we 


have a relation between the concentrations a 


and y, namely 
t O(Y,) OF Yn (12) 


” 


uw ( Tr.) 


In the second optimum problem we have to 
AW, q (% 
N 


v . 57... 

q a | vy 1 Vn Ai " q — *n (Y, A) 
! ! 

since g ts constant we can only maximize the sum- 


maximize QQ, A w,, 


N 


ry) 
and 
mation. Accordingly let 


\ 
Max 2'v,, (y, 
ro} (1 


Ly (a) A) where a, a (13) 


Then the functional equation corresponding to 


(10) ts 
Ay} (14) 


Ly (a) Uy," 


where 


U4 = O(%) (15) 
and 


(16) 


£o (a) 0 


Clearly v,(y, — A) = 0 when v, = 0 or y, = A. 
As v, increases y, decreases until y, 
that a > d(A)) 


A) is a continuous positive function of v, 


A (pro- 


vided, of and since 


Ui" 
between these values it must have a maximum in 


course, 


this range. 

The simple case of a linear equilibrium relation 
r — ay can be solved analytically and shows the 
steps of the process most clearly. In this equation 


(15) becomes 


mM xY 


hence 
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£, (4) 


Max v, (y, — A) 


Max °, 


i 


By straightforward differentiation we find that 


provided a A %, 


£, (a) 


when 


Again 


| 
£, (a) Max Ake 


t 


a 
- A) 
z 


and substituting from (17) and (1S) and differen- 


tiating we have 


2s (a) 


when 


x, and a, (xAa*) (vt) 


is obtamed for the optimal 
two-stage is particularly By 
the the tank ts 


using the optimal one-stage policy with respect 


The way in which v, 
poles imiportant, 
principle of optimality second 
to its in putu (xAa* ’ then 


by (19) vr, 


yee by (21) r 
1/2 , is 
(ar, A) 7 (x* a A) x vy. 


Thus optimal policy is to divide the wash water 
equally between the two stages. 


It is quite easy to see that in general 


(N 1) (x* 


Ly (a) 


when 


x, ands, (xAa’>)' wre (233) 
Further the optimal poleyv is always to distribute 
Thus 
whe re 


total 


the wash water equally between all stages 
the total wash is a fraction 
Wy N [(x* a ays st x}. 


wash water W, h. 


Vy oq 
Thus if 


bq and solving for 


tin 


Vy 
A the appropriate value is 


a(x x) 
vt 


the 


(24) 


Thus the solution § to first = problem 


and N. RR. 


AMt Nbeson 


f / N41 . 
(a, 0) qa id )2 J 
. HK y | Ny x 
Nety 
N« ' 
Nq 
To study the case of a very large number of 
let N + x. Then 


staves we may 


A x) CXD (h xq) Vth) 


f. (a, by) qa (1 (1 h xq) exp (h x4) | (27) 


Another approach to case of large N can be muck 
by considering a somewhat artificial apparatus in 
the rate 4 
through a water 

>) dz 


>and 2 dz. 


which rallinate tlows with constant 


bed and an amount of wash 


pisses cross-current between the plan Ss 
If v (2) is the 


stream 


concentration 
the 


solute 


in the ratlinate and y (2) ub [x (2) 


concentration in the wash water, then 


(ZS) 


inlet 


profit function ts 


is the and ', the exit concentration, 


“ hve me te 


now 


re 


4 q- 


This holds 


Wir) ws pul equal tow x we recover the previous 


for any equilbbriam relation ¢, if 


results for a linear isotherm. 


$. An EXamMPuLi 


As an example of somewhat greater difliculty 
and interest than the linear isotherm, the equili- 
brium relation shown in Fig. 2 was used in a 
on the Univae Scientitic 


series of calculations 
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well bedevil similar programmes can best be 


appreciated if we show first what ts required in 


the final result. The structure of the programme 


shown in Fig. 3 will then be obvious and its 


particularly sticky points may be examined. 





The final result is required in the form of a 
table such as Table 1 which gives optimal policies 
for the isotherm of Fig. 2 and for N 
stages with A 0-05. These Tables begin with 
the value a 0-03 which is the next tabular value 
above the critical value a d (A) 0-021. Below 
this value of a it is not worth extracting with 


wash water “ costing” A 0-05. If we require 





the solution to the second optimal problem it can 


be obtained as follows ; for brevity we will extract 





only the three stage optimum. If the inlet con- 
m 0-20 the bottom 
line of the section N 3 of the table gives 


centration to the first stage isu 


immediately that the maximum profit function ts 

g, (0-2) — 0-1075 (column 9) with a total wash 
Computer Model 1103 at the Numerical Analysis water ratio V, 1-19 (column 12), and wash 
Center of the University of Minnesota. Taking a water to first tank of 7, 0-647 (column 11). To 
maximum concentration of 0-2 for the inlet con find the wash water to the other two tanks we 
centration ensured that all concentrations wer use the principle of optimality and interpolate 
below this so that there were few problems of om previous sections of the Table. Column 10 


scaling. The problems which did arise and might — shows the exit concentration from the first tank 


Table l. 
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%. Flow chart for computations, 


is 2 0-09 and this is the feed to the subsequent 


1 
pair of tanks which must be using optimal policy 
the section 
O-O54 
of the 


stage sub- 


with respect to this feed. Thus in 
N 0-09 we tind 4, 
and v, the second tank 


2 on the line for a 
0-284. Since 


3-stage processes is the first of the 2 


process this wash water is the v, of the 3-stage 
process, hence v, = 0-284 and x, — 0-054 is the 
feed to the third tank. Interpolating between 
the entries for 0-05 and x, the 
0-260. This last step 
for 


? 1 OG mn 


0 
section N 1 we find v, 


not by subtraction 


Us V; 
tion is 7, 


is strictly necessary 
vy 0-259. 
0-0328 so that 83-6 per cent of the 


It is of interest to note 


The exit concentra- 


Ve 


solute has been removed. 
that if the same total wash water were divided 


equally between the stages 80 per cent would be 


v4 


extracted and if it were all used in a single stage 
75 per cent extraction would be obtained. The 
gain may look comparatively slight but it must 
be remembered actual protit is gy multiplied by 
the flow rate g and the value of the extracted solute 
and this may represent a considerable gain. 

The aim of the computation is thus to build 
up such a table section by section. This is a very 


th section 


straightforward matter for when the » 
1)"", the latter 
may be printed out and the new section stored in 
The the 
now serve to calculate the next section of the 
table the It 


convenient to keep a column for Vy for this 


has been calculated using the (n 


its location. instructions in machine 


without least modification. is 


expedites the solution of the first optimal problem. 
The flow chart of the programme in Fig. 3 shows 
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PROFIT FUNCTION 9, 








b — 3k). In either case the maximum is 
located to within (2 3) (2a K). Continu- 
ing in the same way L iterations will 
locate the maximum to within (2 3) 
(2a K). Choice of L and K 
made in the light of the function’s be- 


must be 
haviour, but the search routine can be 
programmed for arbitrary L and K. 
The linear case was found to be useful 
for the and 
examining its accuracy. The calculation 


checking out programme 


of one section of the complete table from 





which Table 1 is edited took about three 





3 4 5 
NUMBER OF STAGES-N 


Pr fit function vs. number of stages. 


Fic. 4. 


the way in which this is put together. The major 
computational difliculties arise in the interpolation 
and search routines. In the first place the equili- 
brium relationship is notably flat 
r= O01 and 0-16 and interpolation is difficult 


between 
when it is required to find wv corresponding to y 


of about 0-18, For this reason it is well to express 
(14) as a maximum with respect to variation of x, 


£y (@) 


A | 
Max | | | 
Max én (11) (a a,) [ w (7) } 


) 
( ?, 


For N 


two local maxima, a broad plateau in the middle 


- 1 the expression in the brackets has 


of the range and a sharper but slightly higher 
maximum at = the This 
presents difficulties to a simple search routine of 
the Fibonaccian type, but the following technique 
The 


is sought 


upper end, situation 


value of the 
first 


was found to be successful. 


function whose maximum was 
calculated at A equally spaced points of the 
the found by 


comparison, 


maximum value 
By choosing A’ sulliciently 


40) the function was found 


interval and 
direct 
large (in our case K 
to be unimodal in the range consisting of the two 
subintervals about this maximum. In this sub- 
interval, say (b, b + 3k), values are calculated at 
4, =b+kandb ~ 2k. If the value at the first 
point is greater than that at the second the 
2k). If on the 


2k is greater 


maximum is in the range (b, b 
other hand the value at a, b 
than that at b + the maximum is in (b + &k, 


minutes. Once a programme of this type 

has been set up any other problem of the 

same class can be solved by inserting a new 

equilibrium subroutine and the appropriate A. 

Similar calculations were done for A 1 10, 

1 80, 1 320 and 1 1280 and the maximum profit 
functions are shown in Fig. 4 for N 

If the cost of 


providing and maintaining one stage is a fraction 


and for infinitely many stages. 


pw of the value of the extract, as mentioned in 
Section 2, then it is clearly disadvantageous to 
add a further stage when the jump in gy is less 
Thus 0-2, A= 0-05 and 


than yu. with a 





20 


LOG eI/A 











TOTAL SOLVENT TO FEED RATIO, Vy 


Fic. 5. 





R. Anis, D. F. 

m 0-004 we find g, — g, = 0-005 but g, 
0-0026, hence three stages are better than four. 
Finally 
solved by means of Fig. 5, where Ty is plotted 
If the total quantity 


of wash water available is 6 the corresponding 


the first optimum problem can be 
against the logarithm of A. 


value of A for any N can be found ; the calculation 
is then done for this A and the optimum policy 
the 


available this is preferable to interpolating in a 


and profit obtained. Since programme 1s 
set of tables for various A. In practice a mixture 
the wash 


but be 


of the two problems might arise, Lc., 


water might have a certain price A 
available only in total quantity less than 4. If 
b in the plane of Fig. 5 
given N. then 


clearly for this number of stages the maximum 


the point A YY, Ve 
lies to the right of a curve for 
wash water capacity 4 ts not required but only 
the value 
On the other hand if the pomt lies to the left of 
the curve N the total capacity 4 ts needed but the 


Vy lving on the curve N for A rN’. 


optimum policy corre sponds to a value of A on 


the curve at Vy = 6. This value of A will be 


vreater than A’ so the wash water 


‘scarcity value.” 


acquires a“ 


5. Discusston 


The notion of dynamic programming seems to 
provide the chemical engineer with a very valuable 
tool for attacking the problem of optimizing a 
multi-stage process. The example given here has 
deliberately been made as simple as possible to 
bring out the essential method of approach. It 
that the 


generalized quite easily to cross-current extraction 


however, problem can be 


is clear, 
with miscible solvents. 

In this case more state variables will be needed 
to specify the feed. Suppose that the solute B 
is being extracted from a solvent 4 with a second 
solvent C. Let a 
Bi (A B) in the raflinate and extract streams 
respectively and s and S the fraction C (4 — B) 
If a 
emerging from the x“ stage and q is the tlow rate 
of (A B) in the mainstream, u the rate of 
(A B) in the extract is the rate of 
the pure solvent C, balances of A B. B and 


and y denote the fraction 


in these streams. sullix n denotes values 


while % 


C respectively give 


and N. R. 


AMUNDSON 


q n~l q n 
Inn 


Yn 1*n-1 Vn Sn 


If equilibrium is attained at each stage there is 
a set of three relations between the four quantities 
r, & y. S, such as is normally exhibited in a 
Janecke diagram. Thus if « be given s(x) is 
known from the raflinate portion of the equili- 
brium curve and at the other end of the tie line 
y (a) S (7) 


relations, three mass balances and three equili- 


and are known. There are thus six 


bria, to determine the six quantities q,, u 


The 


n* - 


n* n* 


Sas Yn » raphical construction for these is 


nn’ 
well known. 

By analogy with what has gone before we might 
seck to maximize the value of the total amount 
If A is 
as a fraction of the extract value. 
N 
2 (u,, Yn A uw n)s 
1 
A w,) can certainly be calculated 


extracted less the cost of the solvent C, 


again this cost 


we seek to maximize The 


vali of (4, Yn 
for any stage when the state Gunes Sa-te Sans of 


, are chosen. Thus when 


the feed stream and x 


have been chosen optimally 


Let 


‘, tt 


the profit function depends on qo. 1. So. 


\ 


Ly \Yo- To Xo) Max 2 (M,, Yn 
u ! 


” 


then the principle of optimality gives 


Bn-1(%- ry, 4,)| 


Ly (Yo. o- So) Max [u, 4) AW, 
«, 


Here it seems that the “ curse of dimensionality ” 
has again descended for gy a function of three 
variables would have to be calculated and stored. 
However, we observed that except for the first 
stage where s, — 0 the s and x of the feed stream 
the raflinate equilibrium 
fn [ qo: Tq. 8 (19) |. ty 
certainly satisties the same functional relation and 


are always given by 
relation. Thus if fy (qq. 9) 


for the complete N stage process we can finally 


write 


Ly (qq Lg: 0) = Max [ (uy Yy — Aw) + fur (Gp r,)]. 


I 
This has reduced the optimum profit function to 
one of two variables but a further reduction is 
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possible by observing that the mass balance equa- 
tions are homogeneous in the flow rates q. W, U. 
have 


Ls hd ‘sf 
Thus if ww, = 54-1 Ga = KnGe-1 WE 


u (1 


ki) In, and 


Also My, Un ™ {a k,,) Yn Uy Gn 


and all the terms in the square brackets in the 
last equation can be calculated when v, and 
r,-, are known and s,_, is taken to be the equili- 


” 


brium value. Let 


Max (al 


h, (19) 


hy) Av] 


Then ty (do: Xo) = Fo My (%e) 


and we may set up a functional equation 


hy (%%) Max (( hy) uv, Av, ky hy (44) 


and obtain 


Ty (Yo: V9) Yo hy (1o)- 


Thus the problem has been reduced to one of 


calculating and storing a function of a single 
variable. 

Other 
example it might be asked if it were desirable to 
the different 
this case the equilibrium relationships s 


Sr, T) 


generalizations might be made, for 


have stages at temperatures. In 


s(r) 


ete. become = s and there are two 


and T,. <A 


heat balance over the stage will allow the amount 


operating variables at each stage v, 
of heat to be supplied or removed to be calculated 
of this would be included in the 
Another 
incorporate the cost of solvent recovery in a 


and the cost 


profit function. modification might 


more realistic manner. Enough, however, should 


have been said to show the great potential value 


of dynamic programming in chemical engineering. 


NOTATION 


cost of operating a stage with conditions q 


maximum value of profit function in an 
N-stage process with feed stream pp 

number of operating variables at each stage 
number of stages 

protit function 

vector of quantities specifying state of stream 
leaving the n'” tank 

flow rate of first immiscible solvent 


vector of operating variables for nth st 


q (% ty) 
premium on extract of concentration y 
th 


age 


Riy) 
T transformation of stream state by n 
Wn 
«, 7 


stage 


flow rate of wash water to n'" stage 
concentration of solute in raflinate 
concentration of solute in extract 
constant in linear isotherm 

relative cost of wash water 


relative cost of one stage 
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Dynamics of a fluid interface 


Equation of motion for Newtonian surface fluids 


kK. Scriven* 
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The importance of interface rheology to foam stability and to the calming of water 


Abstract 
waves by oil has long been recognized in a general way. Recent research reveals its vital im- 
portance in many free-boundary flows, for example in distillation and liquid extraction. In order 
to understand quantitatively these flows, mathematical descriptions of the dvnamies of an interface 
are needed. 

Herein is derived a completely general formulation of the dynamics of a Newtonian fluid 
interface, that is, one whose rheological behaviour is characterized by its equilibrium interfacial 
tension and Boussinese’s [18] two coefficients of surface viscosity. With consideration of the 
intimate coupling between surface and substrate, the analysis finally vields interfacial boundary 


conditions on two-phase flows 


Résumé— Lo importance de la rhéologie a linterface pour stabiliser les mousses et apaiser les 
oscillations de Peau par de Phuile est reconnue depuis longtemps. Des recherches récentes montrent 
son intérét dans la distillation et lextraction liquide par exemple. Pour l'étude quantitative de 
ces écoulements, il est nécessaire de comprendre la théeorie mathématique de la dynamique de 
Vinterface. 

L’auteur en déduit une form le générale compléte pour linterface d’un fluide newtonien 
dont le comportement rhéologique est caractérisé par sa tension interfaciale déquilibre et ses 
deux coeflicients de viscosité de surface de® Boussinesg. En considérant le couplage intime 
entre la surface et le substrat, lanalyse donne les conditions de la séparation a linterface d’écoul- 


ements a deux phases, 


Zucammenfassung— Die Bedeutung einer Rheologie der Grenzfliche fiir Schaumstabilitit 
und Beruhigung von Wasserwellen durch O1 ist im allgemeinen Sinne lingst bekannt. Neuerdings 
zeigt sich aber auch eine grundsitzliche Bedeutung fiir manche anderen Stromungen mit freier 
Obertliiche, z.B. bei der Destillation und Fliissigextraktion. Zum quantitativen Verstindnis 
dieser Stromungen ist eine mathematische Beschreibung der Grenztlichendynamik unerlisslich. 

In dieser Arbeit wird eine vollstiindige, allgemeine Formulierung der Grenzflichnedynamik 
fiir eine Newtonsche Flissigkeit gegeben, das ist eine solche, deren rheologisches Verhalten durch 
die Grenzflichenspannung beim Gleichgewicht und zwei Koeflizienten der Obertlichenviskositit 
nach Boussines@ur [18] bestimmt ist. Im Hinblick auf die enge Kopplung zwischen Oberfliche 
und Unterlage ergibt die Untersuchung schliesslich die Grenzflichenbedingungen einer Zwei- 


phasenst romung. 


In1 RODULCTION 


between two  homo- 


Tue 


geneous phases contains matter in a distinct 


interfacial region 
physical state : 
exhibits properties different from those in the 
gaseous or liquid or solid states. Accordingly, 
new variables, for example interfacial tension, 
enter into the classical thermodynamic description 


matter in the interfacial state 


of equilibrium systems as soon as interfaces are 
considered explicitly. This is well known; what 
is perhaps not as widely appreciated is that other 
new material properties likewise enter into the 
mechanical description of non-equilibrium systems 
whenever interfaces are involved. Now, in equili- 
brium systems interfaces need not be considered 
explicitly unless the ratio of surface to volume 


*Present address, University of Minnesota, Department of Chemical Engineering, Minneapolis 14, Minnesota. 
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is large, because the contribution of interfacial 
free energy to the total free energy is usually 
small. But the dvnamic behaviour of flow systems 
may be profoundly influenced by interfacial 
effects even though the proportion of matter in 
interfacial regions be extremely small. Further- 
more, motion may even originate in an interface 
in systems that are not in thermal or compositional 
equilibrium. 


When at rest the interface between two fluids 


ordinarily behaves as if it were in a state of 


uniform tension, and it is both possible and con- 
venient to represent the interface mathematically 
as a geometric surface in tension. This representa- 
tion is also appropriate to many flows with free 
boundaries; indeed, it is the basis of the treatment 
of capillarity in classical hydrodynamics. There 
the consideration of equilibrium surface tension 
leads to the conclusion that the normal component 
of fluid stress, or pressure, is discontinuous at a 
curved interface. In classical hydrodynamics it 
is also recognized, in connexion with the calming 
action of oil on water waves, that extension and 
contraction of a surface film produce longitudinal 
the Plateau 

which turn 
the 
interface. 


variations of surface tension 
Marangoni-Gibbs (PMG) effect 
to 
components of thuid stress at the 


give rise discontinuities in tangential 


In contrast with the role of interfacial tension, 
the part that may be played in free-boundary 
flows by other interfacial properties has not been 
widely recognized, These are properties pertaining 
to the resistance of an interface to deformation. 
For example, the departure of the dynamic from 
the equilibrium surface tension in the PMG effect 
can be attributed to a surface dilational elasticity 
or viscosity, whatever its ultimate origin. These 
properties were discovered and much discussed 
qualitatively by scientists nearly one hundred 
years ago. Yet since then they have attracted 
only limited attention and that has come chietly 
from physical chemists interested in two subjects : 
the chemistry of insoluble surface films, which 


often happen to be highly viscous or elastic; and 


the theory of foams, which has long recognized 

in a general way that the PMG effect and surface 

shear viscosity are vital to foam stability. 
Among hydrodynamicists and chemical en- 


oy 


vineers there appears to be a growing awareness 
of the importance of the rheological behaviour 
of interfaces. This evidently stems from several 
recent sources ; (a) studies of circulation inside 
bubbles and drops and renewed interest in 
BoussiNesq’s explanation of their anomalous drag 
coefficients [1]; (b) research on the effects of 
3, 4); 


(c) attempts to fathom the apparent stagnation 


9 


surface-active agents on mass transfer [2, 
of portions of the interface in certain gas absorp- 
tion and vacuum distillation experiments [4, 5, 6]; 
(d) discovery that cellular convection patterns in 
thin liquid layers (Bénard cells) are caused by 


~ 


surface tension variations [7, 8]; (e) evidence that 


drop oscillation and interfacial turbulence in 
liquid liquid extraction are driven by variations 
of interfacial tension [9, 10]; (f) investigation of 
the influence of surface phenomena on interfacial 
area in distillation columns and_ absorption 
towers [11]; (g) analysis of the stability of falling 
liquid films [12]; (h) theoretical analysis of the 
rheological behaviour of emulsions [13], with 
subsequent experimental confirmation [14]; and 
(i) suggestions regarding the action of surface 
transport and surface viscosity in determining the 
stability of foams [15-17]. Thus the time seems 
ripe to begin a mathematical formulation suftic- 
iently general to provide a common description 
of most if not all of the interfaces comprehended 
in this listing. 

to interface 


* two-dimensional molecular world, the dynamics 


DonNan’s reference an as a 
of which is analogous to that of the ordinary 
three-dimensional world of homogeneous phases 
in bulk ”* provides a text for this paper, for here 
we examine the dynamics of substances that may 
be called Newtonian fluids in the interfacial state. 
We assume that the * two-dimensional * molecu- 
lar world can be represented as a geometric 
surface and the material therein as an isotropic 
fluid continuum. We derive in general form the 
equations for a linear dependence of stress on 
rate of strain in the interface. This has already 
been done by Bousstnes@ [18]. 
equations are not in a convenient form because 


However, his 


they are expressed in a special co-ordinate system 


*In his Preface to An Introduction to Surface Chemistry 
by J. W. McBatx, Cambridge University Press, 1926. 





L. E 


coinciding with the principal rate-of-strain axes, 


The simple generalization of Boussinesa’s equa- 


tions to arbitrary co-ordinates in a stationary 


surface has been noted by Ericksen [19] and 
O_proyp [13]. Here we make the further exten 
sion to an interface contained in a time-dependent 
interface undergoing 


surface, that is, to an 


continuous change ith shape and extent a con 
spicuous feature of its actual behaviour. We show 
that the equations of dynamic equilibrium are 
completely analogous to those which apply in 
three dimensions, and 
the 
equations of motion in the interface. 

Thus, to find the two- 
dimensional Navier 
Although 


at first glance this might appear to be a 


with 
the 


then combine these 


stress-rate of strain relation to obtain 


here Is 
the 


Poisson equations of hydrodynamics. 


what is don 


analogue of familiar 
trivial 
problem, it turns out to be complicated by the 
the kor 


tunately the difficulty is easily overcome through 


curvature of two spaces ot mterest. 
the use of well-known theorems from differential 
geometry. 


But the 


sullicent for our purposes; we must take 


two-dimensional analogue is not 
an addi 
tional step without analogue in classical hydro- 
dynamics. Because flow in an interface is perforce 
accompanied by flow in the adjacent bulk phases, 
the final st pis to connect the intrinsic formulation 
of the equations of motion in the interface with 
and accelerations m 


stress resultants, velocities, 


three dimensions. 


Let us 
surtace. 
late 
third dimension except as an abstract mathemat 
Unlike the Flatland 


our have 


imagme ourselves inhabitants of a 


Like the Lines and Poly gons who popu- 
Flatland* we are nat 


\nnorrs aware of a 


notion inhabitants of 
that 


curvature 


ical 


we tind surface may at times 


curvature which changes with time. 


seck to 


material 


As hydrodynamicists and engineers we 


describe mathematically the motion of 


in our world; at the moment we wish to do this 


for a certaim class of thuicds 


With this as 


(two-dimensional) 


which we choose to call Newtonian. 


«\ square ik. A Anpsort), Flatland 1 Romance of 
Many Dimensions, (6th Edition) Dover, New York 1952. 


SCRIVEN 


our aim we begin by considering the rate of strain 
in a continuous fluid medium, and this requires 
a mathematical preliminary. 

We can arbitrarily locate a system of curvilinear 


co-ordinates (v.42) in’ our surface and then 
determine the corresponding surface metric tensor 
a,, such that the distance ds between two neigh- 


bouring points is given by 


(ds)? a,du'du” (2 12: £ 


From the variations of the components of a,, 
with position we can deduce the total, or Gaussian, 
curvature of our surface, and from their variations 
with time we can tind the change of curvature 
with time. Having done this we are able to make 
use of the results obtained by mathematicians 
who have studied the intrinsic geometry of a 
surface, Our progress will be easier if we continue 
tensor shorthand, as 


to use their 


given by 
McConneu {20}. 

Now if ds in equation (1) measures the displace- 
ment of two neighbouring particles of thud at any 
then the 


measured by the rate of change with time of ds, 


instant, rate of strain of the fluid is 


as seen by a Lagrangian observer moving with 
Hence at the 
outset it is convenient to adopt a co-ordinate 
that the lines ua’ 
drawn in the thaid 


with it: for in such a convected, or 


one or the other of the particles, 


svstem wu such 


constant 


are lines and deforming 
continuously 
co-moving, co-ordinate system the uw’’s and the 
time are independent variables, the rate of strain 


bemg measured by [21] 


Dids? Da \y 
du 


du> (2) 
in nt 


where JD) Dt denotes partial differentiation with 
respect to time with the convected co-ordinates 
held constant. From this expression the rate-of- 
strain tensor is defined as 


1 Da, (3) 


Later on we shall require the rate-of-strain tensor 
in a fired co-ordinate system, that is, an unchang- 
ing system whose origin remains at a particular 
Making use of a theorem 
of Otproyp [21] we find that the tensor whose 


point im the surface. 
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convected components are Da, Dt has the fixed 
components 


Da, O49 


3 V* wa 
Dt dt ; 


V" Ais (4) 


A 


where V* is the surface, or intrinsic, velocity 
vector and a comma signifies surface covariant 
differentiation [20], a kind of differentiation so 
defined as to evade the complications of curvilinear 
Thus the 


fixed co-ordinates is* 


co-ordinates. rate-of-strain tensor in 


rue Eat 
CONSERVATION 


STRESS AND ATIONS OF 


MOMENTUM IN A 
SURFACE 

In our simple continuum model the forces acting 
on a (two-dimensional) body are either external 
forces which act on every part of it or else internal 
tractions which act across any line element inside 
the body. The former we denote by the applied 
force vector F*. The latter we characterize by 
a stress tensor 7" because the vector representing 
internal traction depends on not only the location 
hut the of the 
Then if m, is the unit surface vector normal to 


aa 


given 


also orientation line element. 


gives the 


of 


a line element, the vector T” 
that 


m, 


traction on element for a state 
stress 7’. 

From (the two-dimensional form of) Newton's 
second law we have the principle of conservation 
= of thuid 


according to which the 


of momentum for an arbitrary area 
bounded by the curve C, 
acceleration and force components in the direction 
of 1, are related by 

| | vA" tds 


| F* 1, dZ 4 ie m1, ds (6) 
Ja A 
z z ( 

where L, is any constant, parallel surface vector 
field. The line integral around the area may be 
transformed by means of Green's theorem for 
a surface [20]; we then have 


*Equation (5) can also be obtained directly by taking 
the material derivative of equation (1); see Trorsprn 


(22). 


fluid interface 


[ y4*°+F*+ TT” )UdzZ=—0 (7) 
5? 


since L, , = 0. This equation is true for any area 
and for any orientation of L,; therefore the 
integrand must vanish identically at every point 


in the fluid. Hence: 


F* + T” 


(3) 


This is (the two dimensional form of) Cauchy's 
law of motion [22]. 

According to the principle of conservation of 
momentum the angular acceleration and moments 
about the origin are related by 


A” u 


dz 


. 


| | «,, Few dd p «,, T*’ m,u' ds 


** . 


f 


(9) 


Making use of Green’s theorem, equation (38), and 


the identity u’ , = 5’, we find by a similar 


argument that at every point in the fluid 
0 (10) 
Therefore the stress tensor is symmetric. 


RATE 
Suni 


STRAIN 


‘ 
Fi 


STRESS AND oO} IN 


NEWTONIAN ACI IDS 


We call Newtonian those isotropic fluids in 
which stress depends linearly on rate of strain. 
The most general connexion between the stress 
and rate-of-strain tensors consistent with this 
detinition is 


Yin (11) 


where the tensors P*’ and E*”* are independent 
of the rate of strain and completely characterize, 
respectively, the static and the flow behaviour 
of the thiid. 


tions in which the fluid is locally in thermodynanue 


If we restrict our attention to situa 


equilibrium, then the components of P*’ ana 
k*” are fluid properties that are determined by 
its nature and local thermodynamic state, e.g., 
by its local composition, temperature and state 
of strain. 

Because stress and rate of strain are symmetric 


tensors, both P*’ and E*** must be symmetric 
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with respect to interchanging the indices » and A, 
and ke 
to interchange of the indices x and 
definition a Newtonian tluid is isotropic, that ts, 
all of its material properties at any poimt are 
independent of direction at that point: from this 
it follows that P* and E*” are isotropic tensors. 
This fact greatly simplifies equation (11), because 


must also be symmetric with respect 
8. Now by 


the most general isotropic, contravaniant surface 
tensor of second order is ra“’, which is symmetric, 
and the 
surface tensor of fourth order having the required 


most general isotropic, contravariant 


svmmetries is €« a” a’ ela" a a” a”. 


Thus. 


T° 


ela“ a’ a’ a™*)S (12) 


that three the scalars 


completely characterize the fthuaid 


and we see properties 


r. € and « 
behaviour. Equation (12) is the (two-dimensional 
form of the) Newton Cauchy Poisson law [22]. 
If we resolve the rate-of-strain tensor into its iso- 
tropic and deviatory parts, we have an alternative 


form for equation (12): 


Ye pa pa 
7 ra xa” a 


‘ 
z 


(13) 


€ (a a’ a” a’ 


Here the new coeflicient, « € «. multiplies 
the rate of isotropic strain and « alone multiplies 
the rate of shear strain." 

In a fluid at rest the rate of strain vanishes. 
Hence T*’ ra”, i.e. 


uniform stress. This stress is equal to +, and it 


the tluid is in a state of 


is also the equilibrium interfacial tension o 


defined in classical thermodynamics. Therefore 


(14) 


T oO 


and, since + is independent of rate of strain, this 
equality is also true for a fluid in motion. In 
contrast with +, the other scalars in equations 
(12) and (13) are flow properties. In 
equation (13), we call « the coefficient of surface 
dilational viscosity and «€ the coefficient of surface 
Like equilibrium interfacial ten- 


view of 


shear viscosity. 
sion they depend on composition, temperature 


*Note the distinction between ¢, a physical property, 
and ¢,,, a permutation symbol (ef. Notation). 


and density. However, there are many flows in 
which they are nearly constant and for the sake 
of simplicity in what follows we shall assume that 


they are constant. 


FORMULATION FOR 


Fiucips 


INTRINSIC 


NEWTONIAN SURFACH 


We require the equations of motion in terms 
of the fundamental variables of flow —the surface 
metric and velocity, Substituting equations (5), 
(13) and (14) in (8) we tind after a little manipula- 
tion that in fixed co-ordinates motion is governed 


by 
yA” F* ana 
(d 2a)), + «[a™*V” 


uA pa 
(a) a 


» BA 


gv" ». (4 2a) ,] (15) 


Note here and below that differentiation with 
respect to time does not commute with the opera- 
tions of covariant differentiation and raising or 
lowering indices. However. we do have the ident- 


ities 


and these are used in arriving at equation (15). The 
equation may be expressed in more convenient 
form through the use of the following identities, 
in which the total, or Gaussian, curvature A of 


our surface first appears : 


bi 


ai 


vy ” 


ela ” V™ 3 


a’ (a 2a) J 


cae" [ 1” 
kK I” 


ao a” + (« 


(d 2a)}, 
(a*") , (16a) 
ea" | mn 


«(2K V" 


(d 2a)] , 
a. 


ua ’% 
re Ts 


(16b) 


The equation for irrotational surface flow follows 
immediately from the latter. 
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And this is as far as we can profitably go with 
a strictly: intrinsic surface formulation of the 
equations of motion. We now drop our two- 
dimensional aliases and view the problem in three 
dimensions, where we recognize that our surface 
is imbedded in and moves through Euclidean 
space. Firstly, we see that equations (16) for 
motion in the interface resemble closely the 
familiar Navier- Poisson equations but in addition 
contain terms that reflect the curvature of the 
interface and its change of shape with time. 
Secondly, we see that the intrinsic formulation 
relates only to acceleration, force and velocity 
components lying in the interface and reveals 
nothing about the resultant forces normal to the 
interface, for example. 

Finally, we recognize that there is no way of 
further identifying the external forces which 
appear in the intrinsic formulation, and that there 
is no for the 
analogue of the continuity equation without first 


basis seeking two-dimensional 
examining the connexion of the interface to the 
substrates. 

SURFACE 


CONNEXION BETWEEN 


AND SURROUNDINGS 

Let us now consider our surface in its relation 
An arbit- 
rary for 
surrounding space is (2', 2®, 2°) and its metric 


to its three-dimensional surroundings. 


curvilinear co-ordinate system the 


tensor is g,. The equation locating our surface 


in space may thus be written in the parametric 
form 
(17) 


2, 3) 


7 N'(a, u?, t) (7 1, 


of the surface through space is 


its displacement velocity, whose 


The motion 
described by 
components are X' = )X' ot. The total velocity 
of a particle of surface fluid is the sum of its 
velocity relative to the surface 
system and the velocity of that system, or, in 
other words, the sum of its intrinsic velocity 
(expressed as a space vector) and the local 


co-ordinate 


displacement velocity : 


OU = Xi" + X,'; = X'/du (18) 


How is the velocity of a particle of surface 
fluid related to the motion of the adjoining bulk 


fluid interface 


phases? In the first place, there is no slip at a 
liquid fluid interface under ordinary conditions. 
That is, the tangential component of fluid velocity 
is continuous across any surface representing the 
interface. Consequently total surface velocity 
and the velocity of the contiguous bulk fluids 
must have equal tangential components. — It 
follows that the intrinsic velocity of a particle 


of surface fluid is 
X’) 
X’) 


X’) = g, X, (W? 


(19) 


where W’ and W” are the velocities of the contig- 
uous bulk fluids. 
In the second place, matter is conserved at 


the interface. On one hand the rate of change 


of mass in the interface is given by 
Nite 
\- > o 
eo 
and on the other hand by 


] | [p (U" 


. | | yd Jax 


* 
, 


— 


) yl" m, ds W') nm, 


: p(U' —W')njdz 
is the unit vector normal to the surface. 
making 


where a, 
Equating these expressions and again 
find that at every 


use of Green's theorem. we 


point in the surface 
po (U" 
p(t" 


(yd 2a)+(y V WW") mn, 


W')n; (20) 


If we assume for the sake of simplicity that there 
the 
have the 


is no net exchange of material between 


surface and either substrate. we then 


surface continuity equation 
ad 
2a 


and the kinematic conditions 


(yV*) 


an, Wi =a, W 


These relations and equation (19) together describe 
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the connexion between the surface and the sur- 
Next we 
the relation of the surface to the surrounding 


rounding fluid. must examine further 
space, 

The temporal variation of the surface metric 
may be expressed in terms of the displacement 


velocity. Thus, since [20] 


aig = X,' Xf g; (233) 


we obtain, upon differentiating intrinsically with 


respect to time : 


dg = 2g, X,' X’ (24) 


The shape of a surface is characterized by the 
change in orientation of the unit vector normal 


to the surface. 


1 28 , | 
i 2¢ € ijk A : A 


as the foot of that vector moves over the surface. 
According to Weingarten’s formulae. this change 
in orientation is in turn described by the second 


fundamental tensor of the surface [20 


DyNAMICS OF A NEWTONIAN 


Fiuerp INTERFACI 
The the 


follow from application of the momentum principle 


equations of motion of interface 


in three-dimensional space. The external force 
acting on an element of surface may have any 
orientation and is therefore denoted by the space 
vector F'. The acceleration of the element may 
have a normal component and is therefore 
denoted by the space vector A‘. But because the 
surface consists of Newtonian fluid. the internal 
traction acting along the edge of the element is 
the surface vector T*, which we have already 
Here, however, we view it as the 
X'T’’ m,. If there 


is no convective exchange of momentum between 


encountered. 
tangent space vector T' 


surface and substrates, then in the direction of 
I, these three vectors are related by : 


[| yA' dd || Fil ds { Til, ds (25) 


e 
LM ‘ 


— 


where J; is any constant, parallel space vector 


field whose covariant surface derivative vanishes, 


SCRIVEN 


A familiar argument hinging on Green's theorem 
and making use of the above identities leads us 


to conclude that : 


yA FANT. 4 n,b,,T — (26) 


at every point in the surface. Resolution into 
components tangential and normal to the surface 


vields a 
Xi(- ya“ +47) 
n' | yn, A’ nb 


If we take the vector product of this equation 
with the unit normal to the surface we recover 
equation (8), which governs motion in the surface. 
the 


and is 


The total motion of the interface includes 
motion of the surface that contains it, 
governed by equation (26). Substituting therein 
the surface stress tensor (equation 13) and its 
16), the 
temporal variation of the surface metric (equation 


covariant derivative (from equation 
24). and the connexion between the surface and 


substrates (equations 19 and 22), we have: 


_ pw 


) F' A.’ ae . 
Dt 


(x + 6) Xa" (a NY W,.) , 
«(2A XN! a (XY W)) 

Ne (NY Ws) , 

ZN ee b(n’? WL) 
n'|2He 2H(«-~ oa XW, ») 


Zeb, « ee (X/ W .)] (28) 


where //7 is the mean curvature of the surface. 
The significance of the individual terms here is 
perhaps clearer when the equation is displayed 
in vector notation. Thus if we detine a surface 


gradient operator 


where the a* are the reciprocal base vectors of 


; are the base 


the surface co-ordinates and the ¢ 
vectors of the space co-ordinates, then equation 
(28) becomes : 

*This result closely resembles the equations of equilib- 
rium in the general theory of elastic membrane shells, as 
is to be expected. Cf. paragraph 12.1 of [23] or paragraph 
26 of [24). 
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«) V,(V,- W) + « [2K 


n)-V.(n-W)) + n[2Ho 


where 


(n Vin. n): Von 


If surface density or composition varies over the 
surface there is a gradient of interfacial tension 
producing the stress resultant represented by the 
the With 


arises another 


first term on right. finite surface 


viscosities there stress resultant 
owing to the gradient of the surface divergence 
if the dilation of the surface is non-homogeneous, 
to 
there is uniform motion in a surface of finite total 
the or 
non-homogencous normal displacement the 


Altogether the first three terms on the 


and yet another resultant owing shear if 


curvature, rotational motion in surface, 


of 
surface. 
right of equation (29) represent the net tangential 
The last 


normal resultant, 


resultant. fourth and term 


the 


arises in a curved 


stress 


represents stress which 


interface from the action of 


interfacial tension, and. with finite surface 


viscosities. from both dilation and shear of the 
surface. 
the 


acting on an element of the 


To complete the analysis we consider 


force F' 
interface. This force is made up of two parts, 


external 


The first derives from any external force field, 
such as gravitation, acting on the entire system. 
The second is the net traction, normal and tan- 
gential, exerted by the substrates on either side 
of the interface, and is the more important of 


the two. Thus: 


F' = y# (30) 


where the traction exerted by each bulk fluid is 
related to the tensor describing the state of stress 


Tn. 
j 


next to the interface in that thuid by 7' 
These stress tensors for the bulk fluids may be 
obtained from the velocity and pressure fields 


about the interface : 


2) n'W", 


mh etl ( W, , 


n'p- (A 


un; g W,,) (31) 


)=mn 


If we were to invert equations (30) and (31) in 


(W 


Vi(n- V, 
Vin 


nn-W)—n W) 


2H (x + «) V,-W — 2e(n n): VW) (29) 


(28), the latter describing the dynamics of a 
fluid 
in explicit form the dynamic conditions to be 


Newtonian interface, we would then have 
satisfied by the flows of bulk fluid on either side 
of such an interface. 

If we really were surface inhabitants we would 
be mystified by the external force fields we could 
infer from the motion of our Newtonian surface 
fluid, 


coupling of surface flow with flow in the adjacent 


Now, however, we recognize the intimate 
bulk fluids. Because of this coupling, if we are 
to the 


simultaneously solve for the accompanying bulk 


solve for surface motion we must 


flows. In so doing, we regard the equations we 
have obtained above as boundary conditions at 
the interface between two fluids in motion. FEqua- 
tions (19) and (22) together make the kinematic 
condition that there be continuity of velocity 
across the interface, and equations (28) and (30) 
make the dynamic condition relating the difference 
between the fluid tractions on either side of the 
interface to its physical properties and geometric 
characteristics. 

APPLICATIONS 


SomeE SIMPLE 


In order to illustrate how this general formula- 
tion enables us to find boundary conditions on 
a two-phase flow with a Newtonian fluid interface, 
we apply it to those simple flows in which the 
interface perpetually coincides with a_ planar, 
surface three shapes 


cylindrical or spherical 


often encountered in practice. In each case we 
use that set of orthogonal space co-ordinates in 
which the interface coincides with an .x°-surface, 


for then: 
y x3 


u! constant 


6, . 
(0,0, 1); 


In the 
the interface is sulliciently small that, in compari- 


following we assume that the density of 


son with the remaining terms, the inertia term 


in equation (28) may be neglected. Note also 
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that throughout this section all equations are 
written in terms of the physical componens of 
velocity and tractions. 

For a planar interface we have : 


X constant 


(33) 


which give the discontinuities in the tangertial 


components of viscous traction across a flat 


interface, and 
(34) 


which expresses the requirement that there be no 
net that 
When there is no variation 


normal traction on an interface is to 


remain flat. in the 


z-direction, equation (32)simplifies to 


ew, 


oy 


€) (35) 


a result previously reported in an analysis of 
* interfacial turbulence ” [10]. 
For a cylindrical interface we have : 


XR 


constant 


12k 


From equation (28) we obtain ; 


When the 


equations (36) and (38) simplify to : 


there is no variation in z-direction, 


1 do 


R dd 


) 1 d l d WwW 
. R dd | R dd 


r. 


(30) 


a (« 


e)jl aw, 
R R (ip 


Md 


(40) 


a result obtained in slightly different form by 
STERNLING and Barr-Davip in their analysis of 
the stability of parallel viscous flow with a free 
surface [12]. 


For a spherical surface we have : 


X§'=4; = 
R* sin? 0; gs, — 1 


constant 


VR; K-1 


R sin? 0; bis 
H k? 


Similarly we obtain : 
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1d 


RO 


j; 1 res 0). WI) 
| R sin 0 ny) dd | I 
j2W, l d 
of ame. 
lk R sin 0 dd 
j J SW, d(W, sin @) 
\Rsin 8 | d¢ may) 


(« 


Ht (41) 


1 do 


R sin 0 dd 


l d 
€)§ —— - 
R sin 0 dd 


fF, : 


(« 


; | > (W, sin @) 
IRsin@ | 
2W 
ELE 
LR 
. d(W, sin @) 
(R sin 0 Nye) 


dW} | 
Na dd |) 
1d 


R dO 


2a 


(433) 


2e{ 1 | > (W, sin @) 
R \Rsin 0 


dW yy 
AT dd | 


With no the 
constant, and with W, directly proportional to 


variation in é-direction, with o 


sin @, equations (41) and (43) simplify to 


2x YW, 


R? 0 
(45) 


These are the equations derived by Boussinesq 
in his analysis of the effect of surface viscosity 
on the drag coetlicients of spherical bubbles and 
drops [18]. Note that our coeflicient of dilational 
viscosity « is defined to be half as large as his. 
Thus far we have dealt with situations in which 
the interface is constrained in both shape and 
extent. As a final example we consider a spher- 
ically symmetric two-phase flow in which the 
interface is extended while retaining its shape, 


such as occurs in simple bubble growth. In this 
case we have: 
X* = w! , Ato w= ¢ F* = Re) 
ay, = gn R? sin? 0; g33 = 1 
W,= W, (0, 0,1); 
n’'W,— R(t) 
Rsin@; by 
H IR; kK 


by R:; by. 


”> 
-< 


From equation (28) we obtain; : 


Pc 
- » 
R\ R 
The second term on the right of the last equation 


reflects the additional differential re- 


quired to overcome the viscous-like resistance 


pressure 


of the interface to expansion. 


CONCLUDING REMARKS 


The equations for these especially simple cases 


can of course be derived without recourse to 


tensor formalism and the theorems of differential 
geometry. The use of these powerful aids, however 
has led to a set of general equations from which 
one can write down, in more or less routine 
fashion, boundary conditions on flows in which 
the 


shapes, even when those shapes are changing 


interface assumes far more complicated 


with time. The real value of this result becomes 


clear when one considers certain physically 
important flows in which the shape of the inter- 
face is continuously perturbed about a_ planar, 
cylindrical or spherical form. 

The foregoing analysis might be extended to 
the dynamics of an interface composed of elastic 
surface solid or of non-Newtonian surface fluid. 
Indeed, the work of OLDRoyp [13] suggests the 
possibility of extending the formulation to include 
all types of linear visco-elastic interfacial be- 
haviour, many of which have, in fact, been ob- 
various insoluble 


served in experiments with 


surface films. 
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NOTATION 


determinant of the surface metric tensor, a, 
acceleration of the surface thiid 

second fundamental tensor of the surface 
external force acting on the surface fluid 
external force acting on the system 
space metric tensor 

mean curvature of the surface 

tocal curvature of the surface 

unit vector tangent to the surface 

unit vector normal to the surface 

surface rate-of-strain tensor 

time 

traction acting in the surface 

traction exerted by the substrate 

stress tensor in the surface 

stress tensor in the substrate 

surface co-ordinates 

velocity of the surface thuid 

intrinsic velocity of the surface thuid 
velocity of the substrate 


space co-ordinates 
At space co-ordinates of the surface 
e-system for the surface [20] 


e-system for the space [20] 


“ijk 


ML? 
MT! 
wr-2 
wr-3 
Mi1-2 


Lower case Greek indices refer to fixed surface co-ordinates. 


surface density 

coellicient of surface shear viscosity 

coeflicient of surface dilational viscosity 

bulk density 

equilibrium interfacial tension 

Upper case Greek indices refer to convected surface 
co-ordinates. 

Lower case Latin indices refer to space co-ordinates, 

Overdot indicates intrinsic differentiation with respect to 
time. 

Curet refers to the bulk fluid located on the — nt-side of 
the interface. 

(Quantities pertaining to the bulk fluid and not carrying 
the caret refer to the bulk thiuid located on the n'-side 


of the interface.) 
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Two-phase (gas-liquid) flow phenomena—I 


Pressure drop and hold-up for two-phase flow in vertical tubes 
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Abstract 
when the liquid flow is annular. Published data on liquid hold-up and pressure gradient have been 


An attempt has been made to develop a theory for two-phase flow in a vertical tube 


examined in the light of this theory together with original data obtained on air water mixtures 
in a jin. diameter glass tube. Predictions of the liquid film thickness were found to be accurate 


to 15 per cent. The theory also predicted the tlow conditions under which annular flow gave 


way to slug flow. 


the gas liquid interface. 


Résumé 


vertical quand écoulement liquide est annulaire. 


An approach has been made to an empirical correlation of friction factors at 


Lés auteurs ont développé une théeorie de lécoulement de deux phases dans un tube 


lYaprés cette theorie, les résultats publiés 


sur le * hold-up ~ liquide et le gradient de pression ont été considérés de méme que les données 


originales obtenues avec des mélanges eau-air dans un tube de verre de jin. de diamétre ; lépaisseur 


du film liquide prévue a été exacte a 15° 


ement pour lesquelles écoulement annulaire céde a Téecoulement en chapelet. 


La theorie prévoyait aussi les conditions d’écoul- 


Une relation 


empirique est utilisée pour les facteurs de friction a Vinterface gaz-liquide. 


Zusammenfassung 


Rohr entwickelt, wenn die Fliissigkeit 


ringformig stromt. 


Es wurde eine Theorie fiir die Zweiphasenstromung in einem senkrechten 


Veroffentlichte Werte tiber den 


Fliissigkeitsinhalt und den Druckgradienten wurden unter dieser Theorie gepriift zusammen 


mit eigenen Messungen an Luft-Wasser-Gemischen in einem Glasrohr von 12.5 mm Dmr. Die 


Voraussagen der Filmdicke der Fliissigkeit 


waren auf 15° 


venau. Ks wurde auch eine 


empirische Beziehung fiir die Reibungsfaktoren an der Grenzfliche zwischen Gas und Fliissigkeit 


abgeleitet. 


INTRODUCTION 

Exameces of moving liquid tilms bounded on 
one side by a stationary solid surface and on the 
other by a gas or vapour phase are commonly 
met in chemical plant. This paper concerns in 
particular the case of concurrent flow of liquid and 
vapour (or gas) in a vertical tube, the liquid 
flowing as an annular film on the wall of the tube. 
It attempts to elucidate the relation between 
the bulk flow rates of the two phases, the liquid 
film thickness and the pressure gradient along the 
tube. 

Various procedures have been proposed for 
correlating data on two-phase flow inside tubes, 
all of them empirical. Many of these are based 
on the superficial velocity of the two-phase 
mixture [1-7] and so can not differentiate between 
the various modes of two-phase flow (slug flow, 
annular ete.). 


tlow, A different approach was 


initiated by Marrinecoi and his co-workers [8, 
® and 10] who related the pressure gradient for 
two-phase flow to that of the gas and liquid 
flowing separately at the same mass rate. They 
detined two dimensionless parameters X and 4, 
where Y is a function of the pressure gradients 
for the separate phases and @ relates the pressure 
gradient tor the single phases to that for two-phase 
flow, and claimed that a unique relation exists 
between them. The parameter ¢ can, however, 
be shown to vary with the flow rate even when 
NX remains constant (see discussion on paper 
by Lockuarr and Martrinectir [10]. Further- 
more the degree of turbulence in one phase of a 
two-phase mixture is not uniquely defined by the 
Reynolds number, the interface between the 
phases being subject to the influences of gravity 
interfacial tension. 


and In spite of its short- 


comings this approach has been adopted, enlarged 
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and tested by a number of experimenters with 


varving degrees of success [11-16] ete. On the 


whole the accuracy obtained is of the order of 
50 per cent and where annular flow conditions 
exist, when R, is generally less than 0-5, accuracy 
is poor; deviations (calcuated on the 
X? — R,) of up to 100 per cent are not 
uncommon, 
Another 
(17, 18) deals with 
assuming laminar flow in the two phases, and 


basis of 


approach, due mainly to Lai, 


annular flow conditions, 


relation between pressure gradient, 


leads to a 
liquid hold-up, flow rate and tube dimensions, 
Its usefulness is limited since it is not possible 
to predict either the pressure gradient or the 
liquid hold-up from a knowledge of the tlow rate 
and tube dimensions alone, 

The case of liquid films tlowing downwards 
over vertical flat plates has been treated more 
rigorously [19] and a rigorous treatment of the 
problem of annular tlow on similar lines has been 
attempted. This paper attempts to advance the 
understanding of the mechanism of tilm flow and, 
for the case of annular flow in smooth vertical 
tubes, to suggest a revised method of predicting 
data on the pressure 


the liquid hold-up 


viven 


gradient and tlow variables. 


Velocity 
Gstridution 4 


Sheor stress 
Orstribution 





Flow 


Fic, 1, 


ANDERSON and B. G 


MANT70OURANIS 


Tueory 
For the purposes of the theoretical analysis, a 
model of the flow conditions is conceived as in 
Fig. 1. 
made : 


The following significant assumptions are 


The liquid film is of uniform thickness at 


any cross-section, 


The pressure gradient normal to the tube 


AXIS 15 Zero. 


There is no velocity slip either between 
liquid and tube wall, or between liquid and 


gas at the interface. 


The velocity distribution inside the liquid 
film follows the * universal velocity profile ” 


of VON Karman [20]. 


The physical properties p, and py, are 


constant in the liquid film, 


Velocity profile in the film 

The universal velocity profile has been estab- 
lished for single phase flow either bounded or 
unbounded and correlates a dimensionless velocity 
_a dimensionless parameter 


parameter, a*, with y 


involving distance from the fixed solid boundary, 


Vapour 


Tube wol! 


Liquid film adhering to wal! 





Analytical model of annular flow. 





A. Single profile. B. Double profile. 
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pu* y 
ft 2 
The characteristic feature of the correlation ts 
the use of the “ friction velocity.” 


This mmplies that apart from the physical proper- 
ties of the fluid and the geometry of the system, 


only the shear stress at the wall, +,. is needed to 


a 
represent the operating variables such as volu- 
metric flow, stream velocity, and total pressure. 
This velocity profile was originally derived from 
experiments on flow inside tubes, ie. under condi- 
tions where the shear stress varies linearly with 
at the 


wall, and is intended to give the velocity distribu- 


radius from a value zero on the axis to +, 
tion in the boundary layer adjacent to the tube 
wall. It has been argued [19] that it ts also applic- 
able to liquid films, and this is assumed in the 
present discussion. In upward annular flow the 
shear stress at the gas liquid interface can be 
greater than the shear stress at the wall whereas 
in single phase flow the shear stress is a maximum 
at the wall, and so this assumption must involve 
the 
shear stress at the interface greatly exceeds that 
at the wall. 
profile is used here as the best available approx- 


some error. which will be greatest when 


Nevertheless the universal velocity 


imation to the truth and its assumption leads to 
the prediction of liquid film thicknesses in reason- 
able agreement with those experimentally 
determined. 

The universal velocity profile is given in the 
form of three equations which fit the experimental 
curve over three regions for convenience called 
the laminar sublayer, the buffer layer and the 
turbulent layer. The expressions are as follows : 


| 20) 
laminar; for y 


buffer: 5 << y* 
In (y* 5)) 
t+ 5Siny* 


30, u* 5-5 


3-05 


turbulent: for y* > 2-5 In y* 


If uw is the velocity of a fluid passing through an 
elementary area dA then the mass flow 


W p udA 


Now d. 1 


27 (7. — y) dy 


and y (u pu*)y* 


and hence 
: F je ja 
W | 2n utur|r y 
: , - pu* * pu* 
Taking » and p as constant 
W c 


» . 
=-7 Ty fA Ww 


For convenience let 


SS « 
yy ’ 
and Re* vm L 

PL 
1 


then [Wipe = Je (ies 


My 

In the laminar sublayer equation (1) gives 
iW 25 125 

, 2  ~=3Re* 

W, for the buffer layer, is given by 


: a l 
WW, 5 (30 In 6 


Re* (450 In 6 
e 


S75) | 
yy 
and the corresponding equation for the turbulent 

layer is 


W, y, (2-5 lny* + 3-05) 


2-13 
Yi Re* 13) 


(1-25 In y; 
5730 
wae os O68 

Re* 


The total mass tlow is then proportional to 
Ww Wy + Wy + We 
These results are most conveniently presented 
by plotting W> against y;* with Re* as a para- 


(Fig. 2). y; (pu* ‘p) 
y; (Re*/r,), y; can be calculated from these 


meter. Since yj" 
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graphs if Wr and Re* are known, Re* is pro- 
portional to u* and u*® can be calculated pro- 
vided +, is known, and this in turn is closely 
related to the pressure gradient along the tube. 

The parameter Re* may be regarded as a crite- 
rion of stability of annular flow, since for a 
particular We, v,;* (and hence the hold-up) has 
a detinite magnitude when Re* is large but 
becomes indetinite as Re* decreases. Re* is 
involved in the \; W relationship because of 
the finite radius of curvature of the surface ; 
when r, co (i.e, for flat plates) y;° ts uniquely 


related to We. 


Momentum balance within the tube 


Consider now a section dz of a vertical tube, at 
a distance = from the lower end, in which annular 
flow conditions are established. = is taken positive 
in the direction of flow. If the total pressure drop 
over the length dz is dP, the pressure vradient 
is dP, dz. Since the pressure gradient perpen- 
dicular to the tube axis is assumed to be zero the 
same pressure drop Is experienced by both the 
liquid and vapour. By a momentum balance 
over the length dz the net force giving rise to 


acceleration of the Vapour is given by 


B. G. MANTZOURANIS 


dF, dP, Ay, 
py Ayy dz + u, dW, 


dP, dP, dP, u, dW, 
dz dz dz sai Ayy dz 


where ap, dz is the pressure gradient resulting 
from ae interfacial shear r+, 


1.€. 4, Ayy r, Ay, dz 


Ts 
In the case of the liquid film the relation becomes 


dP), dPy Ay 


T. 29 f., dz Ay, dz + u, du, 


u 


which can be written 


— ers. 


‘ 4 ‘ > 
IP, dP pe 


dz dz 
u, dW, 


gp, Ry 


(3) 
Ay dz 


The force dF (btt br®) necessary to accelerate 
a fluid stream of mass rate W (Ib hr) flowing 


through a cross sectional area A is given by 


dk d(Wum) 














400 800 





16 OOO 20.000 


W" (Re /4) 


u Re* plot. 





Two-phase (gas/liquid) flow phenomena 


where u — mean stream velocity (W) Ap) and 
m — a constant depending on the velocity profile. 
If the velocity is uniform, In the 
range 40 < y;* < 1000, 1-1 to 1-01 
calculated on the universal velocity 
profile. 


In the 


then m l. 
m has a value 
basis of the 


case of the liquid phase 


dF, = dW, 


Mp) § 


Wy G, Axr Gy, 
R, PL 


now 


Axi pr = =Ayepr 


a( R, =! G, Ayr) 


G 
my d 1 yr | R, o 


m L- 1y roy* 
Ry py. 


and so 


d(in R, " 


? 2 32 
Thus dPoy m, G,? d (In G, ) 


(4) 
dz PL R,? d R, PL 


The introduction of the logarithmic expression 
facilitates computation and the 
the 


thickness and density. 


contribution of changing mass flow, film 


Similarly, for the vapour stream 


dP, = m,G~ 


‘ 
— ~ wag (™ 7) (5) 


In the absence of mass transfer at the interface 
between the two phases the terms 


iu, \ (dW , dW, 
(a) ( ie ) and (7 i | 


in equations (2) and (3) vanish and equations (4) 
and (5) reduce respectively to, 


dP oy _ my G,? “(In Gy, 
PL 


dz PL R, 2 dz 


dP, my G,? d G, 
dz yk, +1 Pi Ry, 


R L 
and 


Since there has been no report of simultaneous 


measurements of liquid hold-up, pressure drop 
and mass transfer (by evaporation and entrain- 
ment) the data available have been examined 
on the basis of these simplified formulae, and the 


possible influence of entrainment is discussed later. 


comparison of 


WorkK 


Data of both pressure drop and hold-up in 


EXPERIMENTAL 


two-phase annular flow are available from several 
12, 16 and 21}. 

Additional experiments have 
the authors [22 
experiments on falling films inside a tube using 
<e with detergent at 20 °C 
Table 1). 


3 and described in detail elsewhere 


sources [3, 6 
been made by 


one of and these also include 


water and and 


water at 50 C (see The apparatus used 


is shown in bo 
[22]. 
and hold-up in two-phase flow. The 
CD consisted of a 0-427in. id. 


50in, Its lower end C was joined to a 3-way 


It was designed to measure pressure drop 
test section 
glass tube and 
long. 
brass cock and the upper end to an on-off brass 
cock, 
bore opening equal to that of the inside diameter 
of the 


way 


These cocks were designed to give a full 
test section and were connected in such a 


as to give minimum interference with the 
flow. 

Hold-up measurements were made by isolating 
the test section and letting the liquid drain over 
and measuring its level in comparison 
length CD. This was done by turning 
the cocks C and D through 90°. 
Pressure-drop measurements were made using 


and B 


mounted axially at the top and bottom of the 


valve 
with the 
simultaneously 


air-purged stainless steel capillaries A 


test section. The end of each capillary was sealed, 
the air leaving through two small holes drilled in 
the capillary wall a short distance from the end. 

For the 
separator G was removed and the liquid intro- 
The 


capillary B was left in position to measure the 


falling film experiments the liquid—gas 


duced into the tube over a circular wire. 
pressure at C above atmospheric caused by the 
drag of the falling liquid film on the air core. 

By visual observation, steady state conditions 
the the 


section for all runs and no visible disturbance in 


were attained in first 6in. of ‘calming 


the flow could be detected on passing through 
the 3-way cock. Deliberate slight misalignments 
of the isolating cocks such as might accidentally 
occur did not produce any detectable change 
in the pressure drop reading. 

A complete set of original experimental readings 
for the concurrent upward flow of air and water 
with the detailed calculations 


together are 
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Two-phase (gas, liquid) flow phenomena 


available [22]. These original results are also be calculated directly; the procedure is necessarily 
shown in graphical form in Figs. 5. and 6. iterative and is detailed below. 
The results and calculations for the failing film 1. R, is guessed and y; r, calculated from 
experiments are shown in Table 1. this. 
2. If the pressure gradient is known the value 


Calculation of liquid film thickness of +, is obtained from the force balance equations 


A theoretical value for the liquid hold-up eannot and hence Re* is calculated. If the pressure 


Table 1 





theor. 


Water at 20°C. 62-2 1b 2-48 1b he ft 


O-1se ‘ is4 33°7 B14 28-1 

W175 iTS | 5 , 36-7 T67 
O252 ‘ Tu ‘ a 50-0 ROG 
Os HS 1lo4s 833°. . H2-5 100) 
Oso | 1470 ; 2-4 73-7 1TOSo 
459 2: 1s { S36 11538 
Os20 32-2 2210 =| 2: 2 O14 1204 
O-SS85 30°! 2580 ay 2-5 100O 1260 
Ooh ; vow 107 slo 
Ons a5°3 B510 5 25-2 lli-4 1558 


Water at 


33 1b hr ft. 


Olio 
135 
0-180 
O255 
Os20 
0-390 
0-450 
0-505 
0-566 
Oolo 
0-650 


O-055 
o-10o 
0-126 
0-168 
242 
O-300 
O34 
O45 
552 
O-SS4 
Oost 


23-6 
2A°1 
B+ 
5-5 58" 51-8 
S2: 63-6 
2s20 102-0 SO T68 
5400 116-5 . 85" 
4160 120-8 7 oO. 
$830 1450) 07°: 
5500 1570) 
Gow 167-0 
OB, 
BOS 13+: 3. 12-4) 
OSS 28°: 2 20-6 
182 s2 
Moe “Be. 
SOS 2. 
1212 77° 
1608 | 102-2 
2003 128-4 
2400 143-0 
2800 151-3 
3210 | 1642 “ 990 


> oe tS te me OF 


G25 
GST 
769 
O16 
1007 
1120 


1180 
1216 


1261 


1200 
1317 


Lubrol W. at 23.2 I : 2-21 2-26 lb/hr ft. 


a i oo 
Ta] | & 
~ fre we Ve 


1070 
1166 
1197 
| beatae 

222 


1261 


oe Se 





Water flowing inside a glass tube, inside diameter = 0-427in. 


The quantities multiplied by 10-8 are in Ib, hr? ft?. 
The quantities multiplied by 10-8 are in Tb, hr® ft, 
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gradient is not known, the friction factor +r, py u,* 


must be given ; this will permit the calculation 
~ 2S before. 


3. The value of W 
$. From the y W 


of +, and of +, 
Is calculated 

Re* plot (Fig. 2) the 
theoretical vaiue of y;* is read and compared with 
the one guessed l(y, Ve) Re*}. 

5. If these do not agree, the procedure is re- 
peated with another value of AR, corresponding 
to y, 7, obtained from y,* theory Re*. Generally 
with the second trial value of R, the two values 


of y;" differ by less than 5 per cent, 


DISCUSSION OF EXPERIMENTAL RESULTS 


The * double” velocity profile 


All the experimental results available from 
the quoted sources have been analysed. In Fig. 4 
the experimental y,* obtained is compared with 
the predicted value for all the above data with 


- 


the exception of the data presented in Figs. 5 


J 


Theoretical 





and 6, which are discussed later. For complete 
fall on the full 
the 30 per cent deviation lines 


agreement, all points should 
diagonal curve ; 
are also shown for reference. 

It is clear that for low values of y,* the predic- 
100 


tion is consistently high whereas for y; 


there is a 15 per cent agreement. This observa- 
tion prompts the following argument. 
Suppose that the interfacial shear 7, Is not 
provided by a gas but by a concentric solid rod 
of appropriate diameter moving at the necessary 
velocity upwards relative to the stationary tube. 
Then it is necessary to consider a velocity dis- 
tribution in the liquid based on the distance 
from the two solid boundaries and meeting halt 
wav between the two surfaces (Fig. 1, curve B). 


. the area under 


If this is so, then for a given y, 


the curve u* against y*, Le. W ut dy*, will 


iY increased ; CONVOCTS ly to accommodate a viven 


W (effectively a liquid flow rate) a smaller yy 





Experimental y* 


Pickrell 
Govier et al, 
Ishin et al. 


Author : 
Oo Karewatani 
Meklwee 





100 


Experimental y* 


Fig. 4. 
locus for single profile 


Comparison of experiment with theory. 
locus for double profile 


Key 

: air water 
air water detergent 
falling film water 20 ¢ 
falling film water 20° (¢ 


falling film water 50 (¢ 


detergent 


116 





Two-phase (ygas/liquid) flow phenomena 





af 


Woter rote 
W000 Ib/r 


Liquid hold-up 











Air rote, W, Ib/nr 


i. 5. Adiabatic air-water flow, volumetric liquid hold-up. Glass tube dia. — 0-427in. Upward co-current flow. 





Woter rote 


J a 


Pressure gradient (d,/dz), (Ibft/nr?)/ft° ft x10 











Air rate, MW, lb/hr 


Upward co-current flow. 


Fic. 6, Adiabatic air-water flow, pressure gradient. Glass tube dia. = 0-427in. 
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is required. The percentage increase in W for 


a given y* ts found to be a maximum of approx- 
imately 12 per cent in the range of y,° 20 — 80, 
The 


limiting Re* values is shown in Fig. 2. On this 


resulting WW y, relationship for large 
basis the locus for complete agreement in Fig. 4 
is represented by the dotted line, and the bulk 
of the experimental points deviate from it to the 
Another result of the 


extent of 15 per cent. 


‘ 


application of this “double” profile is that in 
addition to having the interfacial velocity common 
to both phases the velocity gradients approach 
equality, a condition compatible with the notion 
of equal vorticity of both phases at their boundary. 


It can be argued that a free liquid surface under 


shear stress is not damped to the same extent as 
one in contact with a solid boundary, and the 
velocity profile will lie between the two extremes 
postulated. It is impossible merely from hold-up 
measurements to test this plausible consideration 
because the experimental scatter is greater than 
the resulting change in W*. 

In examining the more accurate and extensive 
results given in Figs. 5 and 6 on a similar basis, 
(Fig. 7) the 
may be made 


following significant observations 


The poimts are arranged in seven groups 
corresponding to the seven liquid tlow-rates 


investigated. 





Locus for doutte profile 


fmeoretical 


y 


+ Locus for single profile 


10% 








5 


y expermento! 


ia. 7 


Letters and figures refer to points plotted in Pigs, 5 and 6, 
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The points of these groups show a scatter 
which is highest at 20 per cent for the 


lowest y;* and lowest at 


6 per cent for 
the highest y;* investigated ; this is as ex- 
pected, although the order of scatter is 


somewhat higher than that anticipated. 


These 


within a band which crosses both the loci 


groups of points can be enclosed 


for single and double profile. 


The scatter in each group is not haphazard 
but shows a definite trend as the air flow 
rate is increased. 


From these observations it may be concluded 
that the apparent agreement or disagreement with 
either profile and the scatter within any one group 
are not altogether due to chance or to experimental 
error, and so further possible explanations must 


be considered. 


Liquid entrainment 


Another possible explanation of the deviation 
of the experimental and predicted values of y; 
is based on considerations of liquid entrained 
in the vapour core. No detailed information is 
available giving liquid entrainment as a function 
of the flow variables. That the extent of entrain- 
ment can be considerable is reported by ALvrs 
[4] who determined this quantity in a limited 
number of experiments involving the concurrent 
flow of air and water or oil in a lin. i.d. horizontal 
tube. Further experimental work on this problem 
has now been completed and will be published at 


a later date. 


The nature of the solid boundary surface 


The relations so far used in connection with the 
universal velocity profile apply only when the 


—?— 
rhe 


solid surface is hydraulically smooth [23]. 


experimental data available are also contined to 


smooth tubes. When the solid surface is rough 


it is found that the expression for the turbulent 
laver, y* > 30, becomes: 2-5lny* + C where 
the constant C is a function of the roughness and 
has a value always less than 5.5. It follows that 
for any given y* the value of u* is smaller when 
when it is smooth. 


the surface is rough than 


Conversely, for a given W* the experimental 
value of R, and hence of y;* will be larger than 
that predicted for smooth surfaces. 

Unless the value of C is known for a given sur- 
face, no analytical treatment is possible. However 
it is clear that the discrepancy between theory 
and experiment already discussed cannot be 
attributed to rough tubes as it occurs in the op- 


posite direction. 


The magnitude of flow rates 


The theory developed deals only with annular 
flow, Le. when no liquid bridges the tube. Hence 
it is essential to be able to determine the cases 
when annular flow breaks down to slug flow. 

Detailed calculations at a constant liquid flow 
reveal that + 


im 


continuously decreases as the air 
until it 


two large numbers. 


small 
The 


resulting value of Re* is rapidly approaching its 


flow decreases becomes very 


as the difference of 
critical value (Fig. 2) and the outlined theory 
breaks down. Near this region it is considered 
that liquid bridges the tube; this is supported by 
This 


air flows well below the point of minimum pressure 


visual observation. condition occurs at 
gradient (on a plot of pressure gradient versus 
gas flow at constant liquid flow). The fact that the 
calculated and theoretical values of y;* agree 
below this point suggests that the minimum value 
does not indicate any change in flow mechanism 
as has previously been supposed. The minimum 
can, in fact, be explained simply by the change 
in the relative importance of the term gp, R, in 
equation (3). This is further supported by the 


absence of any discontinuity in the friction 
factor Reynolds number relationship (Fig. 11). 
In the region of very low liquid flow rates 
combined with very low air rates, after the critical 
value of Re* has been passed, the hydrostatic 
component of the pressure drop as deduced from 
the hold-up measurements exceeds the total 
The 


stress at the wall has negative values despite the 


pressure drop measured, resulting shear 
net positive flow of liquid; this indicates that the 
liquid moves downwards near the wall but up- 
wards away from it, a phenomenon sometimes 


called © slip-back ” in slug flow. 
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Effect of ripples 
It is a qualitative experimental fact that the 
but is 


liquid vapour interface is not smooth 


disturbed by small waves. If these waves are of 
such a nature that the individual particles of the 
surface of the 


circular paths [24], then in addition to the mean 


liquid describe approximately 
surface velocity u;, the particles on the crest will 
have a forward velocity equal to the backward 
velocity of those in the troughs. For tlow over 
a flat surface, the net effect of the ripples on the 
flow rate will be zero. 

In annular flow conditions and generally in 
liquid films, the amplitude of the ripples is of the 
same order of magnitude as the liquid film thick- 
ness [25] and the motion of the ripples may 
the one described. 


deviate considerably from 
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Any deviation would be such as to attribute to 
the crests a forward velocity greater than the 
backward velocity of the troughs because of the 
closer proximity of the latter to the solid boundary. 
The net effect is that for a given mean arithmetic 
the wall, the 


relation to 


film thickness and shear stress at 
vreater the surface disturbance in 
the film thickness, the greater the liquid flow rate 
that can be accommodated. 

As this deviation is expected to be greatest in 
the region of small tlow rates, it is conceivable 
that such an effect might account for what has 
been explained already on the basis of a double 
velocity profile. However, the deviation of the 
experimental y,° from the predicted value is the 
same for the runs with air-water as with those 


when a surface active agent was added to the 
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liquid, although in the latter case the friction 
factors are greater and the visual wave pattern 
markedly different. Furthermore, the same devia- 
tion trend of the experimental y,° is observed 
irrespective of the gas liquid tube combination 
or the flow variables. Hence it may be said that 
the ripples have no effect on the mean arithmetic 
film thickness, at least as far as it is experimentally 
detectable. 

However, the ripples have a very marked effect 


on the vapour phase by modifying to a large 


degree the friction factor +, pypuy, This was 


shown by Laimp [17] who measured pressure 
drops for air flowing in a rubber hose on the 
surface of which waves were artificially induced. 
(See also [11, 25. 26. 27). 

For a complete prediction of the characteristics 
of annular flow, the friction factor for the vapour 
phase must be correlated with known flow vari- 
ables. By making the analogy with rough tubes, 
it is expected that the friction factor is a unique 
function of the Reynolds number of the vapour 
reckoned on the dimensions of the Vapour core 
and of the roughness of the interfacial surface. 
Thus a rational approach to this would be to 


attempt to relate the friction factor to the Rey- 
nolds number and the wave characteristics and 
also to relate the wave characteristics to the flow 
variables and the physical properties of the liquid. 
The realization of such a scheme presents great 
difficulties mainly because even for steady state 
conditions the ripples exhibit a spectrum of 
frequencies and amplitudes. What follows suggests 
an approach to a purely empirical correlation 
between the friction factor for the vapour core 
and known flow parameters. 

The degree of frictional interaction between 
the vapour in the core and the liquid in the 
conveniently expressed as the 
the 
conventional lines, the establishment of a rela- 


annulus can be 


shear stress at interface, 7,;. Following 


tionship between the friction factor 7; p, u® and 
the Reynolds number of the vapour, Re, 
these dimensionless 


py UD, wy is sought. In 


parameters the velocity u might be taken as 
relative to the tube, to the interface. or to the 
velocity of the ripples. 

Since little is known about the ripples the 
choice lies between the first two possibilities, and 
be estimated, both bases were tried. 


as uu; can 
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Cross plot of Fig. 8. 





Since the correlation ts not greatly improved by 


using «,, and as considerably more labour ts 
involved in using it, the simpler basis of correla- 
tion is finally selected. 

To allow for the effect of the liquid flow rate, 


the parameter 


is used, It is calculated from We yo a 
the (VW uf, Re*) plot, 
used in preference to Wy or W 


viven in 
This parameter is 
bre LUNG 

1. When Fig. 9 is constructed from a cross plot 
of Fig. 8 keeping Re, constant, there is an approx- 
imately linear relation between In (+, py u,*) and 
Q*:; the lines for constant Re, the 


slope. This is very useful for interpolation and 


have same 


also for limited extrapolation. From a cross plot 
a tinal set of curves is produced in Fig. 10 for 
integral values of *. 
2. Sets of curves —such as those of Fig. 10 
calculated from the published sources mentioned 
(See Fig. 11), a curve for a 


are such that when 


o| 


lA 


Friction factor 
° 
2) 
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particular value of @* from one worker's results 
is fitted, by a shift of the Re, axis, to the corres- 
ponding curve derived from a different source all 
the corresponding curves of these two groups 
coincide, 

If Q* fails to act as an efficient empirical 
parameter, the same must be expected of others, 
such as the amount of liquid flowing, the ratio 
of liquid to vapour flows (i.c.. in effect, the 
Martinelli parameter), or the local liquid hold-up, 
because there is no a priori reason why these 
should express any property of the ripples. 

The positioning of any one set of curves relative 
to the Reynolds number axis appears to be a 
function of tube size, surface tension, and perhaps 
liquid viscosity. It appears that the bulk velocity 
f the vapour might be a better parameter than 
the Reynolds number for representing the type 
of ripple induced, since the curves derived for 
air water mixtures at atmospheric pressure are 
spaced along the Re, axis by factors corresponding 
to the tube diameters. When the surface tension 


of the liquid is changed either by the addition of 
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surface active agents or by a change of temp- 
erature, the sets of curves can be made to coincide 
by shifting them along the friction factor axis. 

however, insuflicient 


There is, experimental 


data over a wide range of these variables to 


justify an attempt at correlation along these lines. 


Falling films 


In the case of downward concurrent tlow, the 


theory derived for upward flow is applicable 
except that the term gp, in equation (2) and the 
term gp, R, in equation (3) change sign. 

The data of Bercruin [2] which do not include 
measurements of R, have been analysed on the 
basis of this theory; Le. R, was predicted and 
the friction factor Reynolds number relationship 
When the that its 


bulk velocity relative to the tube is equal to the 


calculated. gas rate is such 
interfacial velocity, there is no frictional inter 
action between liquid and gas and hence the pres 
direction of tlow is 
that the 


Berceuin et al. 


sure gradient in. the ZeTO, 


This does not mean friction factor is 


necessarily small, overlooked 
this point and as a result their values for the 
friction factor decrease at low Reynolds numbers 
of the gas. 


When 
liquid drags the 


there vas, the 


with it, 
the gas travelling downwards relative to the tube 


is no deliberate tlow of 


; stationary gas cor 
but upwards relative to the liquid surface; as a 
result the static pressure increases in the direction 
of liquid flow. Thus, in the force balance equations 
the quantities +,, (dP, dz), go, and gp, R, change 
(N.B. the direction of net liquid flow is 
taken as positive throughout this paper). 


Sign. 


The experimental results obtained are presented 
in Table 1, 
single velocity profile basis. The comparison of 
experimental and theoretical y, 
cluded in Fig. 4. 

The locus for agreement with the 


which includes the calculations on the 
values is in 


velocity protile when +; opposes the flow of liquid 


(as in these experiments) lies below the diagonal 


of Fig. 4, since this profile now predicts a greater 


y;* than the single one. The points for falling 
film conditions in fact lie about this locus at the 
higher values of y,*, when the interfacial shear 


is also high (see Table 1). At lower values of shear 
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the pots lie above the diagonal but it is reason- 
able to that, 


gravity-induced ripples have an overriding effect. 


suppose under these conditions, 
Braver 


[28, 20] to account for the influence of ripples in 


Recently an attempt was made by 


falling tilms on the momentum and mass transfer 
rates by measuring the amplitude and frequency 
of the ripples. A. statistical technique was used 
and satisfactory explanations were reported when 
caleulations were based on the so-called residual 
film thickness, ic. after the amplitude of waves 
was taken into account. His interpretation must 
be viewed with some suspicion because it over- 
simplified the actual phenomenon, As his ex- 
periments were arranged so that practically no 
interfacial shear developed, his findings have no 
bearmg on conditions during annular upward 
two-phase flow, 

For tlow on the outside of vertical tubes it is 
found that the curves in Fig. 3 should be displaced 
to the right of the limiting Re* curve along the 
We axis by the same distance that they are now 
displaced to the left. One direct consequence of 
this is that for the same lquid flowrate a thinner 
lilm is obtained than the corresponding one 
inside the tube of the same diameter, an effect 
enhanced by the smaller negative interfacial shear 


stress generally experienced for unbounded flow, 


CONCLUSIONS 


1 OA 
experimentally determined is developed for pre- 


method involving no constants to be 
dicting the liquid hold-up given the pressure drop 
and flow variables for the case of annular upward 
two-phase flow in smooth vertical tubes. It proves 
successful when applied to a wide range of 
experimental data. 

2. It is shown that the flow patterns fall into 
two categories, namely annular flow where the 
core is occupied by the gas phase, and non- 
annular flow where liquid bridges the tube and 
interrupts the gas core. The theory holds only 
for annular flow and predicts the point where this 
type of flow breaks down, 

3%. The presence of ripples on the liquid surface 
is shown not to affect significantly the prediction 
of the mean arithmetical film thickness. 


t. The prediction of the pressure gradient, 
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given the flow variables, depends on the successful Ag = shear area per unit length of tube 
D — diameter 


correlation of the friction factor between the liquid 
KF force required for acceleration 


and gas phases. It is argued that this is a function 
: : y acceleration due to gravity 
of the Reynolds number of the gas core and the } superficial mass flow density 


effective hydraulic roughness of the liquid surface ; ratio of momentum weighted velocity to bulk 


hence a detailed study of shape and size of ripples velocity 
is necessary. However, some indications are given — 
: ‘ : . volumetric flow rate 
as to a possible line of approach for an empirical 
radius 
correlation, hold-up, fraction of cross-section occupied by 
5. Various configurations involving the tlow one phase 


of liquid films past smooth solid surfaces are velocity 
examined as possible applications of the theory, friction velocity at the wall = (7, p,)'/? 
and the success of the theory in the case of falling ee — 
2 ‘ : ; = distance from tube wall 
films inside vertical tubes is experimentally distance measured along the tube in the direction 
established. of liquid flow 


6. The results obtained can be used to refute dynamic viscosity 
a statement which is sometimes made [30, 31), density 
; , surface tension 
and which if true would render the proposed theory 
si 7 . T shear stress 
unnecessary. This is that * the shear stress at the 
4 ‘ : Dimensionless groups 
wall due to a liquid film is the same as that for a Re 


Reynolds number 
full bore pipe with the same average velocity, py, ¢* Y wy, ~ distance parameter 


On this basis, the predicted shear stress at the u u* — velocity parameter 


wall for two runs selected at random is smaller Wy, 2rr, wy (= Re, 4) = liquid flow parameter 
Re* PL u*® Ve M7. 
Q Wey; 


Subscripts 


than the experimental one by a factor of 2 and 


% respectively, 
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Abstract 

are described. The method is capable of giving a true gas temperature with an estimated ac curacy 
of | 2 per cent of the absolute temperature, without the necessity for applying correction factors. 
The temperature measurement can be made nearly independent of any mechanism of heat transfer 
to or from the detecting elements, and therefore may be more accurate than other devices in cases 
where gross errors can be caused by heat transfer phenomena. Such instances may arise where 


very small flows are involved, or low gas velocities exist 


Résumé. auteur décrit les movens de mesurer les témperatures de gaz par utilisation de 


deébimétres seuls on en série, Cette méthode peut donner une température exacte du gaz avec 


une precision estimeé a 2°, de la température absolue, sans qu'il soit nécessaire d’appliquer 


des facteurs de correction. Cette mesure de température peut étre faite indépendamment de 
tout mécanisme de transfert de chaleur entre le gaz ct les appareils de mesure et peut done étre 
plus precise quavee dautres dispositifs dans le cas ou des phénoménes de transfert de chaleur 
occasionnent des erreurs importantes, ce qui se produit pour de trés petits débits ou de faibles 


vitesses de waz. 


Zusammenfassung. Der Verfasser beschreibt Verfahren zur Messung der Gastemperatur mit 


Hilfe von StrOmungsmessern, die einzeln oder in Serie angewendet werden. Man erhilt eine 


wahre Gastemperatur mit einem geschitzten Fehler von + 2 Prozent der absoluten Temperatur 


ohne Verwendung von Worrekturfaktoren. Die Temperaturmessung kann nahezu vollig 


unabhingig vom Mechanismus der Wiarmetibertragung zu oder von dem Anzeigeclement ausge- 


fihrt werden und diirfte daher genauer sein als bei anderen Geraiten, bei denen ein vrosser Fehler 


durch die Wirmeiibertragung entstehen kann. Beispiele dieser Art sind sehr kleine StrOmungs- 


mengen oder niedrige Geschwindigkeiten. 


Tur problem of measuring the true temperature general way of carrying out such determinations 


of a flowing gas frequently becomes rather of gas temperature based on the use of flowmeters. 


complicated, particularly when the total gas flow is | For example, if the mass flow rate is known, then 


very small, or when the gas velocity is low. In a flowmeter can be used to determine a kinematic 


such cases, various temperature sensing elements viscosity, which is, in turn, a unique function 


may be used, and an attempt made to obtain For this case, 


tri 


of temperature for a given fluid. 


gas temperatures by employing radiation it is obviously desirable to select a type of flow 


shielding, high velocity jets, or combinations of — clement with a which is sensitive to 


both the these 
precautions it is still possible to have extremely 


response 


around element. Even with temperature changes, for example, a capillary 


meter operating in the region of Poiseuille flow. 


large errors in the measured gas temperature. The necessity for operating such a device at a 


Caleulated corrections may be applied but these controlled constant mass rate of flow can be 


are also subject to a very considerable error. avoided by using two such devices in series. This 


A measurement of gas temperature which can 


be made much less dependent on the mechanism 


of heat transfer to a sensing clement would be of 


value in these special cases. This would be true 
even though such a method gave only moderately 


accurate results. This article describes one 


arrangement may take two forms. Identical flow 
clements could be used, one of which has gas at the 
unknown temperature passing through it, with the 
being maintained at a known 
Alternatively, 
different flow elements could be used in series, 


second clement 


controlled — temperature. two 
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providing the flow behaviour through each shows 
a different dependence on temperature. Then, in 
the 


direct 


pore sSSUT« drops 


case, by 
the che Vicous, a 
this and the 


flowing gas can be established. 


either comparing 


ACTOSS two relationship 


between ratio temp rature of the 


devices are described in the following 


Although it is not a 


Such 


sections. precise means of 


measuring temperatures, it ts not subject to gross 
errors, and will always give the true gas tempera 


ture within the accuracy of tts calibration. 


THEORY 
Capillary Flowmeter 
A capillary operating in the region of streamline 
low where Potscuille’s law ts obeved will closely 
follow the equation 
aE Doh px 


wah 


W (1) 


where W ois the mass rate of flow, is the head 


across the capillary expressed in terms of height 


Db and L ure 


respectively, 


of manometer thud of density py, 
the capillary diameter and length 
and « and p, the viscosity and density of the 
flowing thaid. 

For a capillary of tixed dimensions, and for a 
specific manometer fluid, with the units of pressure 


drop as specified, this equation can be written ; 


where A, is a constant for the particular capillary. 
\ single measurement of pressure drop and tlow 
rate at a known temperature (¢.g. room tempera 
ture) is suflicient to determine A,, providing the 
capillary dimensions are not seriously affected by 
If the 


measurement of 


of flow Is 
the 


temperature changes. mass rate 
known, for example, by 
flowing gas before it is heated or cooled, then 


equation (2) can be written 
fA A, h. (33) 
py W 


The term yp, is the kinematic viscosity of the 
gas, which is a well known and fairly sensitive 
function of temperature for most gases. It is 
therefore possible to determine the kinematic 


viscosity, and thus the temperature of the flowing 


was, if the mass rate of flow and pressure drop 
across the capillary are known. 

Two capillaries in series. Uf two capillaries are 
used nm series, one at the unknown ius tempera- 
ture 7. and the 
after bringing it to 
temperature, T. then the 


holds since the mass flow ts identical in cach case : 


other using the same gas tlow 


some known controlled 


following equation 


(4) 


- 1 Py er) 


fa T, 


where the subseripts 1 and 2 refer to the individual 
capillaries. Inasmuch as all conditions are known 


and constant at the reference temperature, T 


equation (4) can be written ; 


If A.! is determined at the known temperature, 


then, for capillaries following Poiseuille’s law, 


the temperature at any other point can be deter 


mined from the ratio of drops across 


pressure 
the capillaries and the known kinematic viscosity- 
temperature re lationship of the gas. 

If the capillaries do not behave ideally, then 
the above equation can be modified by treating 
it as an empirical relationship, and writing, 


” 


‘ 
(' } (6) 
My’ T 


If one assumes a simple exponential relationship 


between kinematic viscosity and temperature, 


equation (6) can be modified to : 
(7) 


In this determimations of the ratio of 


pressure drops at two or more temperatures are 


cause, 


required to evaluate the two constants x and m, 
The assumption of a simple exponential relation- 
for the should 


with 


dependence 


fairly 


ship temperature 


apply good over wide 


accuracy 
temperature ranges. 

Orifice and capillary in series. The necessity 
for bringing the gas stream to a known tempera- 
ture can be avoided by the use of two dissimilar 
For example, an orifice or 


elements in. series. 


nozzle has a different temperature dependence 
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with respect to pressure drop than does a capillary. 
One form of the equation for meters of the oritice 


type may be written [1] as: 


. ! 
22 py Pw ne) 
l gp 


W e¥S, 


(5) 


where ¢ is the flow coeflicient, Vo is the so-called 
‘expansion factor,” S, is the area of the orifice, 
and A, is the pressure drop across the orifice in 
units of head of manometer fluid of density pay. 

If the orifice to be used is restricted to small 
diameters, so that it conforms to the case of a 
small nozzle in a large diameter pipe, and if the 
pressure drop across it is small relative to the 


absolute pressure, then equation (8) becomes 
W Ky. (py ho) (9) 


The constant, Ay, in this equation applies to a 
specific nozzle and manometric fluid. 

If a capillary and a small nozzle are placed in 
series, then equations (2) and (9) can be combined 


to give 


Koy (py ho) (10) 


or 4 (11) 
hob? 

According to this equation, the ratio of the 
pressure heads across the two flow elements is a 
unique funetion of temperature for a gas of a 
given composition, It can be seen, also, that this 
combination gives less sensitivity to temperature 
The constant, A. 
temperature 


than does a capillary alone. 
little 
over fairly wide ranges. 
follow their 
sufliciently well, A 
single measurement at a known temperature. 
Again, if the behave 


ideally, one can write equation (11) in an empirical 


variation with 


If both capillary and 


should) show 


nozzle respective flow equations 


could be determined by a 


flow elements do not 


form. 


(12) 


In this case. a number of measurements are 
required to effect a calibration. The exponent, a, 
ean be determined the 


h. to hy at a single temperature for varying flows. 


from 


relationship of 


If measurements of h. ho" are made at two or 
more temperatures, then the relationship between 
temperature and pressure drops can be established, 

The empirical forms of the equations given in 
equations (7) or (12) are probably the most useful 
for accurate work, since factors such as the loca- 
tion of pressure taps, entrance and exit losses, etc., 
could cause minor deviations from the standard 


capillary or orifice equations. 


EXPERIMENTAL TEMPERATURE 


MEASUREMENTS 


Single capillary. This application has been briefly 


described in an earlier article [2]. Very small steady 
flows (100-500 cm? min) of nitrogen and hydrogen were 
measured at room temperature in calibrated rotameter 
type gas flowmeters (Matheson Co. Flowmeters Types 
202, 205 and 204). 
Soft of 1 OD. 
inside an air bath oven equipped with an air circulating 
The held 


within 1 FF. and was measured with both thermometers 


The gases were then passed through 


Sin. copper tubing which was coiled 


fan. oven temperature could be constant 


and shielded thermocouples, 

As the heated gas left the oven it passed though a 
capillary made of a 3in, length of pure silica capillary 
Imm bore. Pressure drops 


tubing of approximately 


were measured on an ordinary water-filled manometer. 


The arrangement for heating the gases was more than 
adequate to bring the gas within | or 2 °F of the measured 


oven temperature, 
Fig. 1 gives a comparison of oven temperature 
with the temperature calculated from the capillary 


pressure drop and the known rate of flow. The 
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Oven temperature, 7 °F 


1. Temperature measurement using a single 


flowmeter at known flow rates. 
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Fic. 2. 
thent 


nozzle in series 


determined at 
It can be seen that the 


oven temperature and the calculated gas tempera 


capillary constant was room 


temperature in this case, 
ture within 5 F over the 
70 #F 
the gas tlow rate and holding it at a steady value. 


A device 


in order to test the applicability of equation (12) 


agree range of 


500 °F. The major error lies in measuring 


( apillary and nozzle. was constructed 
Such a temperature measuring device should be 
as small and compact as possible, and the capillary 
and nozzle should follow their respective equations 
These 


requirement that the capillary must operate at a 


reasonably well. conditions impose the 


orities 
still 


desires easily measurable pressure drops of small 


much lower Revnolds number than the 


at the same mass tlow rate. However, on 


A. woter 


Fic. 3. 
nozzle, h 


Relationship between pressure drop across 


9 and pressure drop across capillary, h.. for 


manometric thermometer 


Manometric thermometer 


of gas temperatures using capillary 











for measure- 


and 


but similar magnitude across each element. 


together with the minimum physical dimensions 
which are convenient, One method of accomplish- 
ing these ends ws shown in the sketch of the device 
given in Fig. 2.) A disk of porous sintered glass 


(Corning pyrex “ medium ~ grade) was used as a 


capillary, thus giving a reasonable pressure drop 


and flow rate while maintaining a close approxima- 


tion to Poiseuille How (i.e. flow through the dise 
followed Darey's law closely h. The nozzle was a 
standard orifice from a small natural gas pilot 
burner, having an opening 0-O18 in. in diameter. 
Before ente ring the capillary dise, the ius passed 


through another fritted glass piece which acted 


‘ 


Fig. 4. Calibration of manometric thermometer. 





Flowmeters for the measurement of gas temperatures 


as a filter. Flow rates used were of the order of 
150 450 ml min which gave pressure drops in 
the range of 1-10 in. of water for cach element. 
Ordinary manometers were used with no special 
taken to 


readings. It is readily apparent that the physical 


precautions obtain highly accurate 


dimensions of a suitable device can be varied 
quite widely by matching porous dises of varying 
sizes and or permeabilities with any suitable 
nozzle. 

In Fig. 3 is shown a logarithmic plot of orifice 
head versus capillary head for two temperatures, 
using air as a test gas. The lines are parallel 
within the experimental accuracy, with a slope 
equal to the exponent “a” in equation (12). 
The of the 


with temperature is shown in Fig. 4. 


ratio of flowmeter heads 
This chart 


variation 
was prepared by plotting logarithmically the 
pressure drop across the capillary agaist tem- 
perature for a constant orifice head. The slope 
of the straight line obtained is a little less than 
the 


measuring system was not completely isothermal 


expected from theory, probably because 
at the higher temperatures. However, the validity 


of equation (12) is apparent for this case. 


Discussion or Errors 


Without special care, readings of manometers 
should be possible to 
This that the 
ratio as given by equation (12), should be obtained 


1-2 per cent accuracy. 


means necessary pressure drop 


with a maximum error of | 3 per cent. If special 
precautions are taken, a maximum error of 
1 per cent could be realized. Allowing for some 
inherent error in describing gas properties as 
simple exponential functions of temperature, the 
maximum error for temperature determinations 
should lie in the range of 2.3 per cent of the 
absolute value, depending on the care used, and 
on the properties of the fluid in question. This 

that the temperature of a 
1,000 “R 


20 F by a single reading. 


would mean gas 
determined to 


Several 


flowing at could be 
within 
readings would probably reduce this error by a 
considerable extent. 


The flowmeter clements are best applied by 


being placed directly in the line carrying the gas 


flow, either as a part of the line, as a by-pass 
arrangement, or as an insert in the duct or line 
with venting to a lower pressure. 

One field of 
temperature measurements of this type would 


veneral application in which 
be useful is that of high temperature gas reactions, 


either catalytic or non-catalytic in nature. 
Temperatures, flow quantities and pressures in 
these studies are dictated by kinetic (and some- 
times diffusional) considerations, and combina- 
tions of conditions frequently arise where the 
measurement of .rue gas temperatures becomes 
a major problem. For example, a temperature 
the kind here 


has been used successfully to give gas tempera- 


measuring device of described 
tures in diffusional studies using low flow rates 
of gases of known composition at temperatures 
up to 300 C, when other temperature measuring 
instruments had shown very large errors. 

Other types of flow elements for example, 
critical orifices, will give the same results, and 
when these are properly combined will yield the 
same kind of calibration plot. While the measure- 
ment of temperature by this technique is probably 
useful only in specialized situations, it can afford 
a positive means of giving gas temperatures 
with good accuracy in such cases, and without 
the 


of any kind, 


necessity for applying correction factors 
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NOTATION 
constant exponents 
diameter 
gravitational constants 
pressure drop across capillary 
pressure drop across nozzle 
numerical constants 
capillary length 
absolute temperature 
mass rate of flow 

= gas density 
manometer fluid density 
gas viscosity 


numerical constants 
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Abstract—The adhesion of a spherical particle to a liquid-air interface has been investigated 
experimentally, and equations derived for the force and work required to detach the particle 
from the interface as functions of the solid and liquid densities, the surface tension and the angle 
of contact. The critical centrifugal force required to detach and submerge siliconed glass beads 
(radii 0-02-0-1 em) from air liquid interfaces was in good agreement with theory. The liquids 
used were water, and 10 per cent and 25 per cent ethanol in water. Similar measurements with 
aggregates of approximately equal sized beads showed that the adhesive strength of the aggregate 
decreased as the number of beads increased, reaching a limiting value of about 50 per cent of the 


adhesive strength of the individual beads, when ten beads were present in the aggregate. 


The bearing of the results on the process of attachment and the adhesion of mineral particles 


to air bubbles in froth flotation cells is discussed 


Résumé— Des recherches expérimentales sur Vadheéesion de particules sphériques a interface 
air liquide ont permis d’établir des équations relatives a la force et au travail nécessaires pour 
détacher la particule en fonction 

Des densités du solide et du liquide 

De la tension de surface 

De Tangle de contact 


La foree centrifuge critique nécessaire pour détacher et submerger les perles de verre aux 


silicones (rayons 0,02-0,1 em) de Vinterface air-liquide est en bon accord avee la théorie, Les 


liquides utilisés sont Peau, et 10°, et 25°, déthanol dans Peau, Des mesures analogues avec des 


agrégats de perles de dimensions a peu prés égales ont montré que la force adhésive de lagrégat 
diminue quand le nombre de perles augmente, atteignant une valeur limite denviron 50°, de 


la force adhésive des perles individuelles quand il vy a 10 perles par agrégat. 


L/ auteur discute lincidence des résultats sur Padhesion de particules minérales aux bulles air 


dans les cuves de flottation par écumage. 


Zusammenfassung— Die Adhiision eines kugeligen Teilchens an einer Grenzfliiche zwischen 
Flissigkeit und Luft wurde experimentell untersucht. Die Kraft und die Arbeit zur Abrtennung 
des Teilchens von der Grenztliche wurden in Gleichungen dargestellt als Funktion der Dichten 
des Feststoffs und der Fliissigkeit, der Obertlichenspannung und des Kontaktwinkels. Die 
kritische Zentrifugalkraft, die zum Abtrennen und Versenken von silikonisierten Glasperlen 
(Radien 0,02 bis 0,1 em) an der Fliissigkeit-Luft-Grenztfliche gebraucht wurde, war in guter 
Cbereinstimmung mit der Theorie. Die verwendeten Fliissigkeiten waren Wasser und 10- und 
25% ige Athanol-Wasser-Mischungen. Ahnliche Messungen mit Aggregaten von angeniihert 
yleicher Teilchengrésse zeigten, dass die Adhiisionskraft des Aggregates mit steigender Teilchenzahl 
absinkt bus zu einem Grenzwert von etwa 50°, der Adhiisions kraft des Kinzelteilchens, wenn 


das Aggregat aus 10 Teilchen besteht. 


Die Uhertragung der Ergebnisse auf das Haften von Mineralteilechen an Luftblasen in 
Schaumflotationszellen wird diskutiert. 





INTRODUCTION 
Avrnovuen the surface chemistry of processes for 
the concentration of minerals by froth flotation 
has been the subject of numerous experimental 
investigations, few quantitative studies have 
been made of the mechanics of the attachment of 
the solid particles to the air bubbles [1], (2). [3]. 
The ease with which the particles can be attached 
adhesion, 


to the bubble and the strength of the 


are however of prime importance when 
effectiveness of a flotation process ts considered, 
When a hydrophobic particle is attached to a 
flat air liquid interface, it is in a condition of 
equilibrium between the gravitational force on 
the one hand, and the sum of the surface tension 
and buoyancy forces on the other. Its potential 
energy, as a function of its height with respect 
to the liquid surface is a minimum. To submerge 
the particle a force must be applied which must 
The magnitude 
the 


density and shape of the particle, as well as the 


exceed a certain critical value. 


of this critical force is a function of SIZE, 
surface tension of the liquid and its angle of 
the 


principles apply toa particle attached to a bubbk 


contact with particle surface Similar 
in a flotation cell, except that the foree due to 
the must considered, 
The the the 


strength of adhesion of the particle to the bubble, 


internal pressure also bx 


magnitude of critical foree. i.e. 
is then a measure of the flotability of the particle 
and it could be used to assess the effectiveness 
of the 


energy 


flotation process. A similar potential 


barrier is involved when a particle ts 
moved from the liquid through the interface into 
the gas, and the magnitude of this also is important 
when considering the process of attachment of 
particles to bubbles in flotation cells 

The object of the work reported here was to 
demonstrate the existence of the critical force, 
to compare its value with theory, and to demon 
how it varies with the above mentioned 


For 


reasons the measurements were carried out using 


strate 


variables. theoretical and experimental 


spherical particles. Silicone-treated glass beads 
were used, and the critical force determined by 
subjecting the particle and supporting liquid to a 
centrifugal field which was increased steadily 


until the surface tension and buovaney could no 


. Nuri 


longer support the particle and the latter sank 


under the centrifugal force. 


THUEORY 
Consider the system shown in Fig. 1 with a 
spherical particle, centre A, supported by surface 
at a flat interface ECBD. Let the whole 
horizontal plane about a 


tension 
system rotate ina 
point O, such that the radius of rotation, R, is 
large compared to the radius of the particle, r. 
The centrifugal force, fF. acting on the sphere is 
given by: 

, ar pla 


3 


(1) 


where p, is the density of the sphere, and a is the 


centrifugal acceleration. 


O 


. 


1 Hy chropobicd vic surface 


tension al an 


sphere supported by 


vir liquid interface. 


The total surface tension force FP. acting around 
the perimeter, can be expressed by : 
F (> 


2 cos 8) (2 rr cos x) 


Ie r 


Jnr ycos x. cos(m7 — 0 — a) (2) 


where y is the surface tension of the liquid, @ is 
the angle of contact between the liquid and the 
surface of the sphere and « is the angle between 
the line joining the point of contact of the liquid 
at the sphere surface and the centre of the 
sphere, and the normal to the direction of the 


centrifugal force. It may be noted here that as 
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Froth flotation : 


x is increased the surface tension component, 
y cos x is also increased, but the perimeter BC 
over which it operates is decreased. By differen- 
tiation of equation (2) it may be shown that the 
force F 


s 


total surface tension goes through a 


when = « 0) 2 and 
0), 2). 

F,. 
from the mass of liquid which would be contained 
in the cylinder BCFG (Fig. 1) together with the 


spherical portion below BC, 


maximum value (7 
4 %e rv cos? l(a 
oman Jarry cos [( 7 
The buoyancy force, can be evaluated 


Hence : 


’ t 
F, 


» Pp, a h a r* cos® « 


py a har cos* x 


7 
r(l 
3 


sin x)? (2 sin «) (4) 
where p, is the density of the liquid, A is the 
height of the evlinder BCFG, and x is the height 
of the spherical segment within the evlinder, 
1 r(l 

An 
derived by arguments similar to those used to 
If the 


effect of curvature of the line of contact of the 


sin %). 


approximate expression for # can be 


estimate the height of a sessile drop [4]. 


interface around the sphere is neglected, then 
it may be shown that 


ap, h* 
Me y (1 — sin 8) 
” 


or h (2y, p a)' *[1 sin (7 0 x) |" ° (@ 


Hence : 


I "y 


py, a \(27 p, a) * wr? cos® a 
{1 


7 
3 


0 x) |"? " 


sin (9 


ra sin x)? (2 + sin x)' (5) 
| 

It follows from equations (2) and (5) that for 
any given system there exists a value of « for 
which the sum of the surface tension and buoyancy 
forces is a maximum. If therefore the centrifugal 
acceleration is slowly increased, the sphere floats 
at the interface as long as 
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F, 


but when the centrifugal acceleration is so large 
that 

F. F, (6) 
flotation is no longer possible and the sphere 
must sink. The critical centrifugal acceleration 


for the sphere is thus determined by the relation 


J 
3 


mr? pla 2rrycos ~cos(7 — 0 x) 


Py a)? zr? cos? x 


le 
p, a \? y 


. ° t , 
sin (7 x) |' * 4 a7 


sin «)* (2 sin n) 


sin x)? (2 sin x) 


a’? r (2 yp,)'? 
0 —a)]'* 


x) 0 


cos? x|1 sin (7 


2 y cos % COS (7 0 (7) 
Exact analytical solution of equation (7) is not 
feasible, but it may be solved by trial and error 
using test values of x and known values of @ and 


y to determine the maximum value of a‘? r. 


The potential energy of the sphere supported 
in the interface and subjected to the gravitational 
force Fe, can be derived by considering the net 
force Fy which acts on it when its centre is at 
a distance H below the level interface : 


Fy Fe 


pi é 


Fy 


es ( 
r= (p 
I 3 Pe 


cos ~ cos (a — 0 x) 


cos? x | 1 sin (7 0 


g(l sin «)? (2 sin «) (8) 


rsina +h (9) 


H 


Fig. 2 


radius 0-lem, density 2-0, at the surface of a 


shows Fy as a function of x, for a sphere, 
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liquid, density 10, and = surface tension 


71 5adyn em which has an angle of contact of 


“4 It shows the 


FY which is the maximum force the system can 


with the sphere. critical force 


sustain in the upward direction without rupture 
of the 
sphere, at a 


and 
the 


for sub 


contact between liequanel interface 


critical value of 2 ». and 


critical downward force Fy necessary 


mergence of the particle, when x 


Fig. 2. Net force, Fy (dyn) ona spol re moan mter- 


face. as a function of the angle. x. between the level 


interface and the line joining the centre of the sphere 
andl the contact with the 


pout of mtertace, 


The full curve in Fig. 3(b) shows PF, asa function 
of HT when the sphere passes through the mterface 
The 


spl rm 


from 

il ’ 
and breaks the 
the 


(negative 


air to licquancl phase, discontinuity at 


is caused when the first contacts 


so that | 


level surface quid rises 


until line of contact makes an angle x 


in sien) determined by 
h rsin %, 


As H decreases further, « tends to zero and FP. 
decreases until, in this special case for which 
0 ow > when H 0. x Oand Fy trig B 
(p, p, 2). Further decrease of H gives positive 
Fy 
This con- 

(Fig. 2). 
and the force attains its maximum value Fy”. 
Further decrease of H the 


break away and the sphere is completely sub- 


values of x and later, negative values of 
(i.e. an upward force on the sphere). 
tinues until x attains its limiting value x 
Causes interface to 
merged. 

When the sphere floats at rest in the interface 


it lies at a depth such that the net force, Fy. 





Ay 


my olved 


Net 
when 
interface. Tl 


forve, Fy (dyn), and work Wo (ergs) 
flat 


of sphere below 


a sphere is moved through a 


distance of centre 


level interface (om) 


The the 
passes through the imterface from air to 
a function of HW by 
integration of Fy. dif ic. from the area under 


is Zero, potential energy WW (ergs) as 
sphere 


liquid, can be derived as 


the curve in Fig. 3(b). The full curve in Fig. 3 (a) 
shows the result obtaimed. An analogous proce- 


dure may be appled to the case of a sphere 
moving through the interface in the opposite 
direction i.e. from liquid to air, and gives the 
dotted curves in Fig. 3 (a) and 8(b). The curves 
the 


critical energy necessary for the detachment of 


shown clearly illustrate critical foree and 


the sphere from the interface in either direction, 


EXPERIMENTAL 


The glass beads used for the investigations were 


treated with silicone oil by the following procedure 


in order to prov ide a water-repellent surface 
The 


were first cleaned in hot chromic acid, washed 


having a known angle of contact. beads 
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Froth flotation The 


repeatedly with conductivity water, drained and 


then dried for at least 48 hr at 150 200 C, 
The hot 


oil which had been previously dried at 100-150 °C 


beads were then submerged in. silicone 


for a similar period of time, allowed to cool and 
stored under the oil until required. Before use a 


bead was transferred toa flat, silicone-treated glass 
plate, and its diameter measured under a travelling 
Only 


from occluded air 


microscope, beads which were round and 
used, The 


which had previously been cleaned in hot chromic 


free bubbles were 


bead was next transferred to the centrifuge 
acid followed by 
bead was then 
until the 
that all 


density 


and the desired test licquancl. The 
washed repeatedly with the test liquid 
tension of the 


surface washings showed 


surplus oil had been removed. The of a 


separate sample of the untreated beds was 


determined with water in a specific gravity bottle. 


8 Cc 


VOSS) /. ZZZ4 











Centrifuge cell and holdet 


The centrifuge cell comprised a evlindrical glass 
test tube, Fig. 
with a rubber stopper. 
hole, D, 
polished perspex, B, 
metal disk, 
central shaft of a variable speed bowl centrifuge. 
of the bead 


could be estimated, without parallax error, from 


t, OS om id. and 5-5 em long, 
It was carried horizontally 
rectangular block of 
dull black 


was mounted on the 


in a bored in a 


mounted on a 


which in turn 


The radius of rotation of the centre 


scales (0-05 in, intervals) which were engraved 


on the top (C) and bottom faces of the perspex 


block. For each measurement the centrifuge 


speed was increased very slowly and the bead 


observed by means of a strobotorch flash synchro- 


repeated washing with water 


adhesion of solid particles to flat interfaces and bubbles 


When the 


and was such that the 


nised manually. speed approached 
bead 
the 


radius of rotation 


the critical value, 


centrifugal forces near 
cell, the 


bead was noted. 


was supported by 
horizontal axis of the 
of the of the 
noted from the strobotorch reading at the 
head from the 
boy liquid. In this 


centre The speed 
was 
was thrown 
of the 


made on 


moment when the 


interface into the 


manner, measurements were several 


beads, 0-02-0-Lem radi, for a variety of radii 


of rotation up to 1-8 em, in three liquids. 
used were, conductivity water 
distilled 


mixed-bed of ion-exchange resin. and 10 per cent 


The liquids 


prepared by passing water through a 


and 25 per cent (by volume) ethanol in water. 
tensions of the 
Nouy 
angles of contact of the liquids were determined 
by the 
plate treated with silicone oil in the 
the beads. 
means of a spherometer, for drops at least 4 em im 


The surface liquids were deter- 


mined using a du torsion balance. The 
sessile drop method [4] on a flat glass 
same way as 


The drop height was measured by 


diameter. 
ResuLrs 


The observed values of the surface tensions, 
densities and angles of contact of the liquids on 


the flat siliconed glass plate were as follows ; 


Pure water ) 71 0.9978 
w° 50-0 dyn cm @ = 80° p, = 0-9813 
25° 30-0dyvn em 0 — 59 p, — 0-9647 


vu 


5dyn cm @ = 85° p, 
_ ethanol 5 
ethanol 5 
The density of the particles was 2-839. 
The critical centrifugal speeds (NU rev min), 


corresponding radii of rotation and bead 


diameters are shown in Table 1 and are illustrated 


in Fig. 5. The 


equation (7). The 


Fig. 5 are those 
dotted 


each full curve are 


full curves in 
predicted by lines 
and below those 


assuming either that / always takes 


above pre- 
dicted by 
its maximum value corresponding to « 90 
or that h Although the 


are in good agreement with equation (7), the true 


0, respectively. results 


theoretical curves may lie higher than those 


shown (and nearer the upper dotted curve) since 
neglect of the particle curvature in equation (4) 
leads to low values of 4. However the true ex- 
values of rNR'? for each 


perimental system 
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Table 1. 


Critical centrifugal speeds for submergence 


Nor 


of beads at an air-liquid interface 





riem) Riem) N (rev min) rNvk 
Beads in water 

OOS | 11K 

O45 


OOo4s 


S255 


O-7T684 


O-485 
0-485 
485 
O-485 
O-485 
O-485 


107s 
1075 
1075 
1Pcl4 
1314 
1-784 
OSS 1-797 
485 
485 
485 
485 


Theory equ 7 


10°, ethanol 
oot 
O97 
142 
485 1-187 
(485 Pz 
485 1 
O4SS 
O-485 
O485 


O-485 


Beads in 
485 
OS) 
O-4S5 


258 
oso 
ass 
Osteo 
ooMw 


O-4S5 ool 


Riem) 


r(em) N (rev min) 


rNyvy Rk 


O-O4SS 
OSS 
O-O4SS 


10355 
uw 


120 


OOS 
0-045 
OMS 
OOM 
oar 


Js0 
re 
V70 
moo 
5 ivi 
5 Io 
lew 


12: 
Oris 


~ 
- 2 


Average 


Theory ; eqn. (7) 


Beads i 
On225 1 
O25 1 tel 
Oo Oss 


Oars 


, ethanol 


O-S065 
OTS 
0 5906 
i al 


OOM mori 


Theory : eqn (7) 





the of the observed 


values, since some of the experimental errors 


probably exceed average 


always give low values, e.g. mechanical vibration 


of the cell, contamination of the interface and 


When 


repeated measurements on the same particle it 


aging of the particle surface. making 
was observed that the values obtained decreased 
slowly, e.g. the data in Table 1 for the 0.0485 em 
radius bead. Similarly the angle of contact on 
the flat glass plate slowly decreased when repeated 
measurements were made on the same surface. 


With all the beads it was clearly observed that, 


at speeds approaching the critical value, the bead 
was supported well below the level interface (as in 
Fig. 1) and that values of A were of the predicted 
order of magnitude. 

Measurements were also made with several 
beads supported simultaneously in an air water 
the tend to 


aggregate together into loosely packed rafts, ore 


interface. Under gravity beads 
bead thick. In one series of measurements up to 
ten beads (0-080 0-084 em diameter) were used. 
The full points in Fig. 6 show the observed 


adhesion strengths (N? R) for the aggregates as 
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hia. 5.) Relation between critical centrifugal speed, 


N, sphere radius, r, and radius of rotation, R 


a function of the number of beads present. The 
open circles show the adhesion strength of the 
vn. bead alone, which was added to the system 
of (n 1) beads to provide the data for the 


aggregate containing » beads. Generally the 
aggregate was detached from the interface as a 
whole. For systems with more than four beads 
the single layer raft became rearranged into an 


ayyregate 


two or three beads thick at speeds 
The 
finally detached by rupture at a neck near the 


aguvregate 


— 


approaching the critical speed. 
uppermost layer of beads. The small pocket 
of air trapped in the voids inside the aggregate 
separated and rose to the surface as the aggregate 
bottom of the cell. The 


moved towards the 


spread of results was partly a consequence of 


different packing arrangements in the aggregate. 
When two or three 


of different 
(0-06 0-09 cm centrifuged — to- 


beads sizes 


diameter) were 
gether, the largest beads were supported deepest 
and became detached and submerged separately 
at slightly less than the speed which was required 
to submerge cach when it alone was present in the 


When different 


were present together the beads formed a loose 


surface. eight beads of SIZES 


. a , 
multilayer aggregate which sometimes detached 
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bia. 6. 
NPR, 


Adhesion strength of agyregates of spheres, 


as a function of the size of the aggregate, nm. 


as a whole, but sometimes only seven beads 
separated in the aggregate leaving one in the 
surface which detached at its own critical speed, 

Some attempts were also made to determine 
the critical centrifugal speed required to detach 
beads from 


glass water into air using a cell 


having a 2mm diameter opening, which was 
completely tilled with water. The results, although 
of the correct order of magnitude, were variable 
and lacked reproducibility. This was attributed 
to the fact that the column of water raised with 
the bead often ruptured and was detached with 


the bead. 


DISCUSSION AND CONCLUSIONS 


The experimental results clearly demonstrate 
that a force or energy barrier must be overcome 
before a spherical particle can be detached from 
a flat interface and submerged. The strength of 
adhesion of a single sphere to the interface is in 
good agreement with theoretical prediction, and 
as expected it increases with surface tension of 
the liquid and with the angle of contact between 
the liquid and the solid. 

For particles having other shapes similar 
principles hold. For some particles, in particular 
those having flat faces, a number of orientations 


at the interface may be possible, each having its 
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particular potential energy barrier and correspond. 
ing adhesion strength. Generally the orientation 
will be such that the line of contact of the interface 
with the particle follows the edges of the particle, 
and it ts necessary to consider the critical depth 
of immersion of the particle such that the line 
of contact can spread around the edge and across 
the Theoretical analysis of this 


plane face. 


problem for a mineral flotation process ts unlikely 


to be practicable. 
The 


importance when the efliciency of a flotation cell 


foregoing principles are of practical 
is considered, Once the particle has been attached 
to an air bubble, strong adhesion is necessary to 
prevent separation during the subsequc nt motion 
of the bubble. 


sudden changes in the rectilinear motion of the 


Inertial forces consequent upon 
bubble can cause separation ; likewise centrifugal 
forces and viscous forces in regions of high shear 
vradients. Only those systems with the strongest 
adhesion can survive and reach the froth laver. 
Direct experimental determination of the critical 
adhesion strength by the above centrifugal force 
technique could therefore provide a useful proce- 
flotation 


and would assist comparison and correlation of 


dure for the evaluation of a process, 


observations on the same flotation process carried 


out in different cells, or in laboratory tlotation 


techniques (e.g. Hallimond tube methods). 
The observations concerning the behaviour of 


aggregates of beads also have an important 
bearing on the performance of a flotation cell. 
Clearly the number of particles attached to the 
bubble and their size distribution are important 
factors when the loss of load from the bubbles 
the flotation 


Aggregates of equi-sized 


and the efliciency of process is 


considered particles 
can be detached from the bubble by centrifugal 
or inertial forces consick rably less (50 per cent) 
than those required to detach the individual 
particles. Mechanical entrainment of small par- 
ticles in voids between larger particles of the 
aggregates can also occur. 

When the foregoing principles are applied toa 
mineral 
bubble. 
force, rr 
the bubble [3] which supplements the 


particle or aggregate attached to a 


some allowance must be made for the 


,. contributed by the internal pressure of 


ravita- 


“u 
~ 


Nor 


tional, inertial or centrifugal foree, and thus 


reduces the adhesion strength. In the case of the 
0-l em radius sphere considered in Figs. 2 and 3, 
F, exceeds Fy” for bubbles smaller than about 
O15 em radius, and adhesion of the sphere to 
This 


size is only slightly greater than the limit set by 


bubbles smaller than this is not possible. 


the lifting capacity (buoyancy) of the bubble, 
For larger bubbles, F, for a 0-3 em 
radius bubble Fo =~ Fy” 2 the effective 


adhesion strength is half that at a flat interface. 


FO eg. 
and 


For smaller spheres the effect is less important 
since Fy’ is proportional to r but F), is propor- 
tional to r® e.g. 
attached to a 03cm radius bubble Fi, = Fy” 4 


i.e. the effective adhesion strength ts 75 per cent of 


for a O05 em. radius sphere 


that at a flat interface. The relative magnitudes 


of f. 


value of the surface tension. 


and Fy” are only slightly affected by the 

It is also of interest to consider the bearing of 
the results on the process of attachment of a 
The 


initial stages of penetration into a flat interface 


particle to a bubble in a flotation cell. 


and the rupture of the surface involves no energy 
Provided H 


always possible for suitable sized spheres having 


or force barrier. r attachment ts 
finite angles of contact and low relative velocities, 
given sullicient time for the liquid tilm to thin 
and rupture. Similar principles hold for the 
attachment of a spherical particle to a bubble, 
except that the work which must be done against 
the internal pressure of the bubble introduces a 
small energy barrier. If the relative velocity ts so 
that the kinetic 

barrier detachment 


large energy exceeds either 


energy for then adhesion 
will not occur, but the particle will either bounce 
back elastically into the liquid, or pass through 
the the the latter 


case the particle will also pass through the 


interface into air side. In 


opposite side of the bubble unless some of its 


kinetic energy is dissipated inside the bubble. 
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Shorter Communications 


An experimental method for determining local mass transfer 
coefficients at a solid surface, and its application to the case of 
a circular cylinder in a transverse air stream 


(Received 6 October 1959) 


INTRODUCTION 


SeverRat authors have described a technique for studying 
interfacial mass transfer phenomena in which the rate of 
sublimation or solution from the surface of a solid into a 
stream of fluid is determined from measurements of the 
dimensional changes of the solid ({1}, [2), [3], [4], [5)). 
This technique promises to be of particular value in 
investigations of the fundamental mass and heat transfer 
processes in two-phase two-component systems, for which 
it offers important advantages over other methods. In 
the first place, solid surfaces are free from the complicated 
effects with the 
surfaces between fluids in relative motion, such as often 


hydrodynamic instability associated 
obscure the interpretation of experimental results for 
gas liquid systems. Secondly, the gauging of local changes 
of dimension makes it possible to determine local mass 
transfer rates at specified points in a solid surface ; whereas 
at fluid interfaces, or at solid surfaces from which the loss 
of material is estimated by weighing or by chemical 
analysis of the fluid phase, only the average rate of mass 
transfer over the whole area can be found. 

In accurate work with this technique the dimensional 
changes which develop in the course of an experiment 
must be kept small enough to leave the flow pattern 
round the solid essentially unchanged. It is therefore 
necessary to use very delicate means of measuring local 
contigurational changes to obtain accurate mass transfer 
data. 


tions is sufficiently sensitive to vield data of high accuracy 


No measuring device so far used in these investiga- 
compared with that attainable by other means. Winpinc 
and Curenery [1] used slip gauges to determine the clearances 
which developed between naphthalene cylinders exposed 
to a transverse air stream and the moulds from which 
they were originally cast. The absolute accuracy of this 
method is limited by the 

O-0025 em, 


available to about 
of 10 
a change 
likely to affect the fluid flow pattern in many cases. The 
[4] and 
Garner and Grarron [5) are apparently similarly insensi- 


gauges 
corresponding to an accuracy 
per cent on a dimensional change of 0-025 em 
photographic techniques used by Frossiimc 
tive. Curistian and Keztos [2] measured the sublimation 
from the surfaces of sharp-edged naphthalene tubes, 
exposed to an axial stream of air, by means of a dial 
indicator capable of registering movements of the measur- 
ing spindle of the order 0-00005 cm. The accuracy of this 


means of measurement is uncertain, however, since the 
pressure of a measuring stylus fine enough to record 
local changes of surface configuration is likely to indent 
the soft surface of the volatile solid. 

In the experiments to be described, the measuring 
means was a fine jet of air impinging normally on the 
surface 
constant pressure reservoir via a restrictor or capillary 


under examination. The air, supplied from a 
construction in the delivery pipe, issued from a nozzle 
placed very close to the surface. Changes in the back 
pressure of the air supply, measured at a point between 
the nozzle and the restrictor, provide a sensitive and 
nearly linear indication of changes in the distance between 
the nozzle and the surface. This technique of pneumatic 
yauging, which allows minute local recessions of the surface 
to be accurately measured without mechanical contact, 
was applied, for trial purposes, to the determination of 
sublimation rates at points round the perimeter of a 
circular evlinder placed in a turbulent transverse stream 
of air. 


EXPERIMENTAL 


The cylinders used were of acenaphthene. approxi- 
mately 1-4 em in diameter and 7-5 em long. Acenaphthene 
was chosen as the volatile solid on the basis of data given by 
Main-Sairut [6], because, (unlike naphthalene), it is not so 
volatile as to sublime appreciably in the gauging jet or in 
the atmosphere of the laboratory in the process of measure- 
ment. Each evlinder was cust from a split wooden mould 
around a steel spindle 1-25 em in diameter and about 25 em. 
long. The spindle was then mounted between centres in an 
ordinary lathe and the end faces of the acenaphthene 
were turned until the cylindrical coating was made just 
long enough to span the working section of the wind 
tunnel used. The evlindrical surface of the acenaphthene 
was also rough-turned in the lathe at this stage. It was 
finally turned true and accurately concentric with the axis 
of the gauging system in the apparatus shown in Fig. 1. 

This apparatus comprised a rigid cast iron stand, having 
uprights A, integral with the base, provided with vees in 
which the steel spindle B could rest horizontally, as shown, 
with one end bearing against the tip of the long horizontal 
lead-screw, C. The turning tool D was mounted on a 
screw-fed kinematically designed slide, KE, having one 
degree of freedom, viz. « sliding motion perpendicular to 
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Fic, 1. 


A Uprights with bearing vees. B- Ground 


kK Kinematic tool slide. F 


mounting. J 


the axis of the work. The cutting edge of the tool was at 
the underside, the work being rotated against it in such a 
direction as to increase the reaction between the spindle 
and the supporting vees. 

The requisite depth of cut having been set by means of 
the feed screw, the spindle and acenaphthene cylinder F 
were rotated by means of a handwheel secured to one end, 
while the lead-serew was slowly rotated by means of its 
knurled knob. 
surface of the acenaphthene cylinder was made to pass 


In a few minutes’ operation the whole 


under the fixed turning tool. 

The same stand and vees were used for the turning and 
gauging, and the tool acted on the surface at a point 
near the gauging nozzle, G. Surfaces generated in this 
way, unlike those turned in a centre lathe, appear almost 
perfectly cylindrical from the standpoint of the measuring 
jet; Le. the 
constant at all positions of the spindle in the supporting 


nozzle surface distance remains nearly 
vees, irrespective of small departures of the acenaphthene 
surface from the straight cylindrical form caused by 
irregularities in the form of the spindle. 
used were silver steel rods of stock size, which are readily 
available ground to close limits of circularity and diametral 
uniformity, and selected for straightness ; 
thene surfaces formed on them could therefore be taken 
to be perfect circular cylinders from the standpoint of 


their behaviour in the wind tunnel. 


The spindles 


the acenaph- 


After turning, the cutting tool was withdrawn and the 


steel 
Acenaphthene evlinder. G 


Circular protractor, I 


Turning and gauging apparatus. 


spindle. (C- Lead screw. 1) Turning tool. 


Gauging nozzle. H — Parallel spring nozzle 


Fixed restrictor. 


from 
The 


was locked in such a position that the mid zone of the 


finished cylindrical surface was freed powdered 


material by wiping with cotton wool. lead-screw 
acenaphthene coating corresponded with the gauging jet. 
The gauging nozzle G was mounted in a flexure device, 
H, essentially a simple parallel spring mechanism furnished 
with a fine-pitched screw bearing against a reducing lever 
integral with one of the stiff parallel Cross-members. 
Using this arrangement the nozzle could readily be set 
suitably close to the surface, to give a high back pressure 
of air. The mid zone was gauged with the jet at 5 intervals 


of rotation of the cylinder, a circular protractor J attached 


to the spindle allowing the angle corresponding to each 


gauge reading to be measured. 

The spindle was then removed from the vees and 
mounted vertically in the wind tunnel in such a manner 
that the evlinder 
horizontal The 
fitted close against the 


acenaphthene spanned the square 
of the 


tunnel walls, 


working section. turned ends 


volatile coating 
and the spindle projected through close fitting holes drilled 
in them. The spindle was set so that the zero division 
of the attached protractor was aligned with the axis of 
the tunnel. 

The air current in the tunnel was induced by a fan driven 
by belting and pulleys from a constant speed motor. 
The air intake section was provided with a paper honey- 
comb and two net gauzes to straighten the flow and 
turbulence ; the tunnel contraction 


promote uniform 
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wud the velocity 
j™ rsy ‘\ 
flat, as 


entrance, ic. at the 


ratio was profile in the 


Tem section Was 
Pitot 


sp tcl 


square, long working 


found to be satisfactorily determined by 
toom the 


mounting holes. A Pitot measurement of the wind speed 


traverses from 


at the axis of the working section and the air temperature 


were taken at the start and finish of each run, and the 


average of these was used in calculating the Revnolds 


number for the 
After 12 


the acenaphthene 


run 
24hr exposure to the transverse air current, 


eviinder and its central spindle wer 
removed and replaced as before in the gauging 
The 
perimeter of the 


but at 


apparatus 


acenaphthene surface was again gauged round the 


yom previously selected. in general at 


5 mtervals, intervals of 1) in where the 


re_htis 
was found to change with angle 
it s« If 


ygauye rapidly 


The jet 


approximately 


pressure 


nozzle was a thin-walled brass tubs 


internal diameter, and 
Negretti and Zambrana 


restrictor, K 


Wiis 


supplied by Messrs tovether with 


the fixed these components are similar 


to those used in the flapper-nozzle assemblies of pmeunnatic 
amd are «ke signed 


control instruments made by this firm 


to operate at an air supply pressure of 180 cm of mercury 
Separate calibration experiments carried out with a bench 
calibrated im ten-thousandths of an inch 
track 


falls nearly 


muacrometer 


imdicated that the pressure of air immediately 


upstream of the nozzle lineariy from TSO te 


approximately 107 cm of mercury as the distance between 


the nozzle and a eviindrical surface of 1-25 em diameter 


is Increased from zero to O4(hRZ19 on The constant 


pressure of the air supply (derived from a compressed 


air eviinder via a reducing valve) and the back pressure 


in the nozzle were measured on ordinary Bourdon gauges 


As these were sensitive to changes of OS om of mereury 


it appears possible (so for as the pneumatic system is 


concerned) to measure changes in dimension over the 


range 0 to O(0R25 em with an absolute accuracy rather 


better than 0-000020 cm, ic. less than 1 per cent of 


the largest measurable change 
In practice, this accuracy was somewhat reduced by 


frictional and hysteresis effects in the Bourdon 


patipves 


used \ further cause of scatter in the observations is 


thought to have been the presence of small and 


this effect might be 


pits 
erevices in the surfaces reduced by 
in creasing the wall thickness of the nozzle at the opening. 
The error involved in removing and replacing the spindles 
in the gauging apparatus appeared to be negligible « ompared 
with these effects 


of relative local rates of mass transfer 


The reproducibility of determinations 
at different points 
round the perimeter of the evlinders was found to be of 
the order 5 per cent 


Resvtirs AND Dtsctssion 


Lacal sublimation rates 


Polar plots of rate of acenaphthene removal vs. angular 
distance from the forward stagnation point are shown in 


Fig. 2. The curves relate to air stream Reynolds numbers 


(based on eviinder diameter) of 10,000, 15.500 and 21,500, 
corresponding to wind speeds in the range 9 21m sec. 
The points on this figure were obtained by converting the 
pressure gauge readings into surface nozzle distances by 
the use of the nearly linear calibration curve obtained with 
the bench micrometer previously referred to, and subtract- 
ing from each of these distances the corresponding distance 
the 
The results, divided by the length of run (and the 
of the 


rates 


obtained before exposure of the evlinder in wind 


tunnel 
solid), vive relative measures of 


Smoothed 


molecular volume 


the sublimation curves were then 
constructed by averaging these relative rates of acenaph- 
thene removal for pairs of corresponding points svmmetri 
eally placed with respect to the wind direction. Curves 
I and II obtained at the 


(TOO) 


were Revnolds number 


The 
75 per cent 


suite 


maximum discrepancy between them is 
similar in genera! form to those obtained 
Rev nolds 
Reynolds number range of these experiments, from: 10.000 


to 21.500, it is seen from Fig. 2 


(1) The 


stagnation 


The curves are 


by other workers at similar numbers, In the 


maximum sublimation rate is at the forward 


pount tines 
that at the reur 


of Powrtiu 


and is approximately 16 
This compares with the ratio 15 


and Geirrrrus [7) and 1-4 of Louriscu 


S|! for Revnolds numbers in the range used here 


The sublimation rate at the forward stagnation 


point is approximately IS times the average over 


the whole surface (obtained by graphical integra 


thon) 
the 


points about SS) from the forward stagnation point 
It is wlant 


The minimum sublimation rate is at pair of 


3040 per cent of the mean value, and 
there is some indication that it increases relative to 


the mean sublimation rate 


with increase of 


Revnolds number 

For the 
spheres in 
YO te Loo 
that the 


flow of water 1v5icm soluble 


the Revnolds 


Garver and Grarron 


puast 


number range from 
5) have shown 
point of minimum mass transfer corres- 
ponds toe the point of flow separation, determined 
near the 
the 


Similarly the 


by visual observation of the flow pattern ; 
the 


alinest 


separation point boundary layer within 


wake appears stagnant 
minima on Fig. 2 occur at approximately the same 
position (85) from the forward stagnation point) 
as the points of flow separation for circular cylinders 
Reynolds 


by Hiewenz, Fact 


in the scann number range determined 


and others [0) 


Rectilinear plots show indications of points of inflection, 
possibly of local maxima, in the sublimation rate curves 
at points between 100° and 120° from the forward stagna- 
tion point. The nature and location of these points is made 
uncertain by the exceptional seatter of the experimental 
data in their neighbourhood. Thev may have to do with 


the influence of the tunnel walls on the flow past the 
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eviinder. it appears that two of Winoinc and Cienry’s 
curves for the same Revnolds number range and a similar 
ratio of tunnel width to evlinder diameter (4-2) show a 


similar inflection, though this is not discussed im |1 


Mass transfer coefficients 


} apour pressure of act naphthe me 


For the calculation of mass transfer coeflicients from 
measurements of the rate of dimensional change due to 
sublimation, data are required for the density and vapour 
pressure of the subliming solid at the temperature of the 
experiment. Ilere the density was taken to be that given 
for solid acenaphthene in International Critical Tables. 
The density of the cast material used may have differed 
from this by a few per cent. The only published data for 
the vapour pressure of solid acenaphthene appear to be 
those given in a paper by Linper (10). These values are 
evidently much too high, however, being of the same order 
as those for naphthiulene, which rough pneumatic gauging 
To check 
this, the over-all mass transfer coeflicient was calculated 
from the data from Run 3 (Curve IID), by graphical in- 


tegration over the whole perimeter of the cylinder, assum- 


experiments indicate to be far more volatile. 


ing Linper’s value for the vapour pressure of acenaph- 
thene ; 
Schmidt number of 2 for the system acenaphthene— air, 


jp was then calculated from this, assuming a 


and an equivalent heat transfer coeflicient was found 
from the assumed equality of jp and jj. The calculated 
result, 2.0 kg cal hrm? ©, 


by well established correlations for heat transfer data from 


is very much below that given 


single cylinders. Thus McApams’ equation 10-2 of [11] 


mass transfer coeflicients at a solid surface 


Fic. 2. Local rates of mass transfer 
from the surface of a circular cylinder 
transverse to a turblent air stream 
Rate of acenaphthene removal as a 
function of angular distance from forward 
stagnation point for 1-4 em diameter 


eviinder at room temperature 


Curves I and II for Re 10.000 
Curve IIT for Re 


for Re 


15.5000 


Curve I\ 21.500 


(Air stream Reynolds numbers 


based on eviinder diameter) 


vives at the Revnolds number of Run 3 a heat transfer co- 
eflicient of approximately 88 kg cal hrm? °C. It is not 
likely thet 
error in the assumed value of the Schmidt number, ie. 


a large part of this discrepancy is due to an 
to an abnormally low value for the difliusivity of ace- 
naphthene ; it thus appears probable that Linper’s value 
for the vapour pressure of acenaphthene is very much toc 


CONCLUSIONS 


dimensional changes due to 


The 


sublimation by 


measurement of 


means of pneumatic gauging has been 
shown to provide a delicate method for measuring local 
mass transfer rates from a volatile solid. The gauging 
apparatus can readily be made sensitive enough to measure 
configurational changes in the range 0-0-0025 em with an 
accuracy of the order 5 per cent; Le. it can be used to 
measure accurately changes small enough to leave the flow 
pattern round the solid essentially unchanged during an 
experiment. 

Measurements of the relative rates of mass transfer at 
different 
eyvlinder placed transverse to a turbulent stream of air 
have been determined by this 


stantially with those of other workers. 


points on the perimeter of an acenaphthene 
means, and agree sub- 
Absolute values of 
local mass transfer coeflicients for this system cannot be 
compared with data from other sources until the vapour 
pressure of solid acenaphthene is accurately known. The 
results of present experiments indicate that the literature 


value is several times too high. 
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Book Reviews 


N.N. Semenov : Some Problems of Chemical Kinetics 
and Reactivity. Vol. 1. Pergamon Press, London 1958. 
305 pp. 50s. $7. 


DurinG the twenty years since Semenov published his 
well-known Chain Reactions the problems of chemical 
kinetics have continued to engage the attention of chemists 
Yet in spite of the vast number of 
that 
SEMENOVY states very definitely in the introduction to the 


almost everywhere. 


papers which have appeared, it is interesting 
present book that, in his view, the main problems concern- 
ing reaction kinetics “ are as yet not completely solved ™ ; 
and this situation he compares unfavourably to the present 
state of knowledge of nuclear reactions and particularly 
chain fission. 

Textbooks for the student may present these problems 
in an oversimplified way, but that is not the aim here. 
It is essentially a treatise for those engaged in research. 
Its subject matter is free radicals and their reactivity, 
and the initiation, propagation and termination of chain 
thus be of extreme interest to the 


reactions. It will 


reaction kineticist but rather less so to the chemical 
engineer (except perhaps for those engaged in combustion 
research). The shortness of this notice is therefore not to 
be taken as detracting from the importance of this volume. 
But indeed a new Semenov book is hardly in need of 


recommendation. 


kK. G. Deneicu 


Grundlagen der Chemischen Technik (Ilerausgeber 
Il. Mohler und ©O. Fuchs) Band IIL- LL. Piariri: 
Werkstoffe der Chemischen’ Technik. Verlag 
Sauerlinder, Aarau und Frankfurt (Main) 1959. S&88s. 
28 DM. 


Das Buch cnthilt in sehr gedriingter Form, klar und 
iibersichtlich angeordnet, eine Fille von Material aus dem 
yesamten Gebiet der Werkstoffe. In einem grossen allge- 
meinen Teil behandelt der Autor in geschickter Auswahl 
der 
Darauf aufgebaut 


die Grundlagenkenntnisse tiber die Eigenschaften 
Werkstoffe und ihrer Oberfliichen. 
kommen die Veriinderungen durch Verarbeitung, auftre- 
tende Fehler, Schiidigungen und ZerstOrungen wihrend 
der Verarbeitung zur Sprache. In einem besonderen Teil 
zusammengefasst wird die Verwendung der einzelnen 
Werkstoffe in der chemischen Technik besprochen. 

Dem Autor ist es gelungen, sein grosses Wissen, verbun- 


den mit longjihriger Erfahrung auf diesem Gebiet, fiir den 
Leser in leichtverstandlicher Art darzustellen. Dass neben 
den vielen Sorten von metallischen Werkstoffen auch den 
Kunststoffen Rechnung wird, soll 
lobend erwithnt werden. Das Buch bietet jedem, der sich 
in das Gebiet der Werkstoffkenntnis will, 
eine rasche und griindliche EKinfiihrung: es sei aber auch 


yebihrend getragen 


einarbeiten 


jedem in der chemischen Industrie tatigen Ingenieur und 
Chemiker, sei er mit Entwicklungsarbeiten oder in der 
Fabrikation beschiiftigt, zur Lektiire bestens empfohlen. 


A. GUYER 


Probleme des Kreuzstrom-Warmeaus- 
Springer-Verlag, Berlin Gottingen /Heidel- 
83S. 10.50 DM. 


Hi. Kou: 
tauschers. 


berg 1959. vii 


Das Bucn behandelt vorwiegend den sogenannten Kreuz 


Gegenstrom Wiirmeaustauscher, Hierunter ist eine 
Anordnung verstanden, bei der einzelne reine Kreuzstrom- 
Wiarmeaustauscher derart hintereinander geschaltet sind, 
dass die HauptstrOmungsrichtung dem Gegenstrom ent- 
spricht. Ein bekanntes Beispiel dafiir ist der Mantel- 


roOhrenaustauscher (“tube and shell exchanger”) mit 
seinen Schikaneblechen. Die Anwendung der Beziehungen 
fiir den reinen Kreuzstrom (Nusselt, 1911) verbietet sich 
dadurch, dass die Temperaturen tiber dem Querschnitt 
nach Durehlaufen einer * Stufe ~ nicht mehr gleichfOrmig 
sind. Der Verfasser berechnet nach von ihm entwickelten 
Verfahren die effektiven mittleren Temperaturdifferenzen 
fiir fiinf verschiedene Anordnungen von Kreuz—Gegen- 
strom Wirmeaustauschern. Darunter ist auch der 
wichtige Fall der spiralig gewickelten Rohre. Bezugs- 
yrésse ist in allen Fallen die effektive mittlere Tempera- 
turdifferenz beim reinen Gegenstrom. Weiterhin behan- 
delt der Verfasser die Einfliisse der Durchmischung in 
den strOémenden Medien und der Wirmeleitung in der 
Rohrwand, die in den meisten praktischen Fiillen zu 
vernachlissigen sind. 

Die Ergebnisse sind in tibersichtlicher Form als Dia- 
yramme wiedergegeben, so dass die Ermittlung der 
effektiven mittleren Temperaturdifferenzen als Funktion 
der Anordnung, des Verhiltnisses der Wasserwerte und 
der Zahl der Stufen leicht mdglich ist. Einige Beispiele 
illustrieren die Unterschiede fiir die verschiedenen Anord- 


nungen. Das Buch ist allen Interessenten zu empfehlen. 


UU. Grieuus. 


Book Recrivep 


J. R. Maver: 


Analysts Pocket Book, Butterworths Scientific Publications, London 1959. 
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NOTICE OF MEETING 





The Protection of Gas Plant and Equipment from Corrosion 


\ joint symposium sponsored by the Institution of Gas Engineers, the Society of Chemical Industry 


Corrosion Group and the College of Advanced Technology, Birmingham, is being held in Birmingham 


on 22nd and 23rd September, 1960. 


Further details may be obtained from Mr. A. 8. Blower, Department of Chemistry, College of 


Advanced Technology, Suffolk Street, Birmingham 1 


Errata 


N. Domprowskt, D. Hasson and D. BE. Warp: Some aspects of liquid flow through fan spray nozzles. Chem. 
Engng. Sci. 1960 12 35. 


Au paragraphe (c) du Résumé la derniére ligne devrait étre: “. . . . la paramétre depaisseur est seulement une 


fonction de (P,)' *  ». 


Die letzte Zeile des Absatzes (c) der Zasammenfassung sollte heissen: ©“ Bei hohen Spritzdrucken ist der Dicken- 


Parameter nur eine Funktion (P,)' * + 


On p. 43 the third line from the bottom of the left hand column should read © (see Plates 22, 23)... .. In Fig. 11, 
the ordinate should read 1 x, cm-'. The caption to Fig. 12 should include : kerosene aleohol O water. 


On p. 44 the eighteenth line from the bottom of the left hand column should read “ at two surface tension levels 
of 30 and 70 dyn em.” The caption to Fig. 14 should read “ The effect of viscosity and pressure on the thickness para- 


meter AK. = Wdyn em O & 7T0dyn em.” 


On p. 46 the second line of the left hand column should read “25-487 dyn em,” 


y4 


PK 10). 


On pp. 45-47 the ordinates in Figs. 15, 16, 17 and 18 should read 


In the caption to Fig. 19, ° 7 aleohol” should read ~ A alcohol.” 





